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Abstract 

to  use.  The  perturbation  procedure  for  the  calculation  of  the  quality  factor  of  the 
composite  resonator  due  to  radiation  into  the  semiconductor  wafer  is  discussed,  aa — * 
no-ted-tha-6  the  perturbation  procedure  enables  calculations  for  the  case  of  rectangu- 
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case  the  calculations  of  the  quality  factor  using  the  perturbation  procedure  are  in  good 
agreement  with  the  results  obtained  from  the  earlier  more  cumbersome  direct  procedure. 
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1.  Introduction 


The  composite  resonator  consists  of  a  uniform  thin  layer  etched  in  a  small 
well-defined  region  of  a  semiconducting  wafer  to  form  a  diaphragm,  upon  which  is 
deposited  a  thin  piezoelectric  film  along  with  the  electrodes  to  form  a  resonant 
region  directly  on  the  wafer.  Under  this  program  the  case  of  the  aluminum- 
nitride  film  on  gallium-arsenide  was  investigated. 

Before  proceeding  with  a  discussion  of  the  work  performed  under  this 
program,  it  is  essential  for  clarity  that  the  meaning  of  the  words  "energy 
trapping"  be  understood.  Since  the  pure  thickness-extensional  resonant  fre¬ 
quencies  are  cutoff  frequencies,  there  is  usually  a  nearby  frequency  range  in 
which  the  transverse  mode  shape  is  evanescent.  There  is  also  a  nearby  frequency 
range  in  which  the  transverse  mode  shape  is  trigonometric.  Consequently,  by  the 
selection  of  the  appropriate  thickness-extensional  overtone  (or  fundamental) 
and/or  the  appropriate  adjustment  of  the  geometry  in  the  electroded  and  unelec- 
troded  regions,  the  transverse  modal  behavior  can  be  made  to  decay  with  distance 
away  from  the  electrodes  in  the  unelectroded  region.  The  resulting  vibration  is 
called  a  trapped  energy  mode,  which  radiates  a  controllably  small  amount  of 
energy  into  the  adjacent  thick  portion  of  the  semiconducting  wafer  and,  hence, 
results  in  the  highest  possible  Q,  albeit  with  many  nearby  spurious  modes  with 
high  Q.  Alternatively,  the  overtone  and/or  geometry  can  be  selected  so  that  the 
mode  does  not  decay  with  distance  away  from  the  electrode  in  the  unelectroded 
region  and  the  resulting  vibration  is  not  a  trapped  energy  mode.  In  this  case 
much  more  energy  is  radiated  into  the  adjacent  thick  portion  of  the  semiconduct¬ 
ing  wafer  and  much  lower  Q’s  result.  Although  there  are  still  many  nearby 
spurious  modes,  they  are  less  troublesome  because  the  Q's  are  lower.  Most 
experimental  work  on  the  composite  resonator  reported  has  been  for  this  latter 


case.  On  the  other  hand,  before  this  program  was  started,  a  detailed  analytical 

treatment  of  the  composite  resonator  for  the  case  when  trapping  is  present 

appeared  in  the  literature*  along  with  a  detailed  discussion  of  when  trapping  is 

and  is  not  present.  Essentially  motivated  by  work  done  under  this  program, 

experimental  results  obtained  at  Westinghouse  when  trapping  is  present  have 

2 

recently  been  reported  in  the  literature  .  However,  the  results  are  for  a 
piezoelectric  thin  film  on  silicon  rather  than  on  gallium-arsenide. 

2.  Discussion  of  Work 

It  has  been  found  that  the  fundamental  essentially  thickness-extensional 
mode  will  not  trap  for  an  aluminum-nitride  film  on  a  gallium-arsenide  diaphragm 
in  the  flat  plate  configuration.  However,  the  fundamental  mode  will  trap  if  the 
gallium-arsenide  diaphragm  is  appropriately  notched  a  small  amount  in  the 
electroded  region,  as  shown  in  Fig.  9  of  Ref.  1.  In  addition,  we  have  found 
that  the  second  essentially  thickness-extensional  mode  will  trap  for  the  same 
film  and  substrate  materials  in  the  flat  plate  configuration  for  a  ratio  of 
film-thickness  to  diaphragm-thickness  larger  than  .69. 

An  analysis  of  the  vibrations  of  a  composite  resonator,  which  is  driven  by 
the  application  of  an  a.c.  voltage  across  strip  electrodes  on  the  major  surfaces 
of  the  film,  has  been  performed.  The  analysis  includes  the  pertinent  waves  in 
the  active  region  of  the  composite  resonator,  as  well  as  all  radiating  waves  in 
the  thick  gallium-arsenide  plate.  The  solution  is  obtained  by  satisfying  the 
differential  equations  for  the  piezoelectric  film  and  semiconductor  as  well  as 
all  boundary  conditions  on  the  major  surfaces  of  the  film  and  semiconductor 
exactly  and  using  the  appropriate  variational  principle  to  satisfy  the  remaining 
conditions  along  the  minor  interfaces  approximately.  The  minor  interfaces 
separate  the  electroded  from  the  unelectroded  regions  of  the  resonator  and  the 


thin  region  of  the  gallium-arsenide  from  the  thick  region.  Both  the  configura¬ 
tion  in  which  the  film  ends  at  the  edges  of  the  electrodes  and  in  which  it 
continues  to  the  edges  of  the  etched  diaphragm  have  been  considered  when 
trapping  is  not  present,  along  with  the  latter  configuration  when  trapping  is 
present.  In  each  instance  the  Q  at  the  resonance  condition  has  been  calculated. 

In  performing  the  aforementioned  calculations  we  have  found  that  it  is 
imperative  that  all  radiating  plate  waves  in  the  thick  region  of  the  gallium- 
arsenide  be  included  in  order  to  achieve  accuracy.  Since  at  a  given  frequency 
the  number  of  radiating  waves  in  a  plate  goes  up  significantly  with  thickness, 
we  have  considered  gallium-arsenide  wafers  no  thicker  than  8  mils  at  a  frequency 
around  132  MHz,  for  which  there  are  30  radiating  plate  waves.  Specifically, 
calculations  have  been  performed  for  thicknesses  ranging  from  1.5  to  8  mils. 

The  1.5  mil  case  was  considered  at  an  early  stage  in  the  calculations  to  check 
the  program  with  as  small  a  number  of  dispersion  curves  as  possible.  All  the 
definitive  calculations  were  for  a  film  thickness  of  7  microns  and  a  diaphragm 
thickness  of  14  microns  and  the  lateral  dimensions  were  adjusted  slightly  to 
maintain  the  same  resonant  frequency  for  computational  convenience.  The  major 
calculations  were  performed  for  wafer  thicknesses  ranging  from  4  mils  to  8  mils 
because  this  is  considered  to  be  within  the  practical  range.  The  calculated  Q 
is  a  very  rapidly  varying  function  of  the  wafer  thickness.  Consequently, 
calculations  had  to  be  performed  for  very  small  increments  in  thickness  in  order 
to  get  all  the  peaks  and  valleys  in  the  interval. 

In  the  absence  of  trapping  in  the  case  in  which  the  film  ends  at  the  edges 
of  the  diaphragm,  the  highest  Q  obtained  is  about  5000  and  the  lowest  is  about 
10  and  there  are  about  10  peaks  for  thicknesses  between  4  mils  and  8  mils.  The 
highest  valley  has  a  Q  of  1000.  The  calculations  were  performed  using  an  incre¬ 
ment  in  thickness  of  1  micron.  In  the  case  in  which  the  film  ends  at  the  edges 
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of  the  electrode,  the  highest  Q  calculated  is  about  15,000  and  the  lowest  Q  is 
about  200.  Calculations  were  performed  when  trapping  was  induced  by  notching 
the  diaphragm  under  the  electrode.  As  expected,  the  Q  due  to  radiation  can  be 
made  as  high  as  we  wish  simply  by  extending  the  lateral  dimensions  of  the  film 
and  diaphragm.  We  have  calculated  Q's  higher  than  200,000  for  quite  reasonable 
dimensions.  In  interpreting  the  foregoing  information  it  should  be  remembered 
that  the  high  Q's  calculated  should  be  higher  than  the  actual  Q's  because  the 
material  Q  and  the  Q  due  to  radiation  into  the  air  are  not  included.  The 
results  discussed  above  mean  that  in  order  to  obtain  reasonably  high  Q  when 
trapping  is  not  present  for  a  given  wafer  thickness,  the  thicknesses  of  the  film 
and  diaphragm  must  be  very  precisely  selected.  The  aforementioned  relative 
stability  and  instability  in  Q  when  trapping  is  and  is  not  present  has  been 
observed  at  Westinghouse  Defense  and  Electronics  Center  but  with  silicon  dia¬ 
phragms  rather  than  gallium-arsenide. 

A  brief  version  of  this  work  has  been  published  as  Ref.  3  and  a  more  com¬ 
plete  version  has  been  published  as  Ref.  4.  Both  Refs.  3  and  4  are  appended  to 
this  report. 

The  type  of  calculation  performed  in  Refs.  3  and  4  is  extremely  cumbersome 
to  perform  and  was  preliminary  to  constructing  a  perturbation  theory  to  calcu¬ 
late  the  Q  due  to  radiation  into  the  bulk  semiconductor,  which  is  much  easier  to 
use.  The  more  cumbersome  direct  calculation  is  required  in  order  to  check  the 
accuracy  of  the  perturbation  calculation. 

A  perturbation  analysis  of  the  Q  due  to  radiation  into  the  semiconductor 

wafer  has  been  performed.  This  analysis  is  considerably  less  cumbersome  to  use 

3  4 

than  the  earlier  direct  variational  treatment  '  and  is  not  restricted  to  the 
case  of  strip  electrodes  and  diaphragms.  In  the  treatment  the  resonant  mode  of 
interest  is  determined  from  the  equation  for  transversely  varying  essentially 


5. 

thickness-extensional  inodes  in  composite  resonators^  and  simple  approximate  but 
very  accurate  conditions  at  the  edge  of  the  diaphragm.  This  resonant  mode  is 
then  used  to  determine  the  radiation  into  the  semiconductor  wafer  by  means  of  a 
variational  approximation  procedure.  Then  the  resonant  mode  and  the  radiation 
field  are  employed  in  a  perturbation  integral  to  calculate  the  Q.  In  this  work 
only  the  configuration  in  which  the  film  continues  to  the  edges  of  the  etched 
diaphragm  is  considered  both  when  trapping  is  and  is  not  present.  However,  in 
the  latter  work,  the  cases  of  rectangular  electrodes  and  diaphragms  are  con¬ 
sidered  . 

3  A 

Of  course,  as  in  the  earlier  work  ’  ,  all  radiating  plate  waves  in  the 

thick  region  of  the  gallium-arsenide  must  be  included  to  achieve  accuracy. 

Calculations  utilizing  the  perturbation  procedure  have  been  performed  for  the 

case  of  strip  electrodes  for  the  same  definitive  geometries  considered  in  the 
3  A 

earlier  work  ’  .  These  geometries  consist  of  wafer  thicknesses  ranging  from  A 

mils  to  8  mils  and  a  film  thickness  of  7  microns  and  diaphragm  thickness  of  1A 

microns.  Although  different  lateral  dimensions  were  considered  in  this  work, 

for  the  strip  case  for  comparison  with  the  earlier  work  a  diaphragm  width  of  500 

microns  was  used  when  trapping  is  not  present  and  600,  when  trapping  is  present. 

The  calculated  results  for  the  strip  case  are  in  good  agreement  with  the  earlier 

3  A 

more  cumbersome  direct  calculations  ’  .  When  trapping  is  not  present  the 


highest  Q's  calculated  are  very  nearly  the  same  as  those  obtained  in  the  earlier 


3,  A 


direct  calculation  *  ,  but  the  lowest  Q's  calculated  by  means  of  the  perturba¬ 


tion  procedure  tend  to  be  nearly  an  order  of  magnitude  higher  than  those  cal- 

.3, A 


culated  by  the  earlier  direct  procedure  ’  .  We  are  not  absolutely  sure  of  the 
reason  for  this  discrepancy,  but  there  are  two  possibilities.  The  perturbation 
procedure  might  be  tending  to  lose  its  accuracy  for  low  Q  because  of  the  in¬ 
creased  radiation  or  the  resonant  frequency  might  not  have  been  sufficiently 


6. 
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3  4 

precisely  determined  by  means  of  the  earlier  direct  procedure  *  for  the  accu- 
|  rate  determination  of  the  lowest  Q  values.  However,  the  fact  that  the  highest 

>  Q's  calculated  by  the  perturbation  procedure  are  slightly  higher  than  those 

t 

i 

calculated  by  the  direct  procedure  seems  to  support  the  latter  possibility.  The 
location  of  the  peaks  and  valleys  of  Q  with  wafer  thickness  determined  by  means 
of  the  perturbation  procedure  is  in  quite  good  agreement  with  those  obtained 
!  from  the  earlier  direct  calculation.  For  the  case  of  square  electrodes  and 

diaphragms,  in  the  absence  of  trapping  the  calculated  Q's  are  roughly  between 
1/2  to  2/3  of  the  values  in  the  strip  case.  This  is  as  expected  because  of  the 
radiation  in  the  two  orthogonal  directions  for  square  diaphragms.  When  trapping 
is  present  the  calculated  Q's  for  both  the  strip  and  square  cases  increase  very 
rapidly  with  the  distance  from  the  edge  of  the  electrode  to  the  edge  of  the 
diaphragm. 

A  brief  version  of  this  work  has  been  published  as  Ref.  5  and  a  more 
complete  version  will  be  published  as  Ref.  6.  Both  Ref.  5  and  a  preprint  of 
Ref.  6  are  appended  to  this  report. 
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I,  Introduction 

The  composite  resonacor  consists  of  a  uniform  thin 
layer  etched  in  a  small  well-defined  region  or  a 
semiconducting  wafer  to  form  a  diaphragm,  upon 
which  is  deposited  a  thin  piezoelectric  film  along 
with  the  electrodes  to  form  a  resonant  region 
directly  on  the  wafer.  In  this  work  calculations 
are  performed  for  the  particular  case  of  the 
aluminum-nitride  film  on  gallium-arsenide. 

Before  proceeding  with  a  discussion  of  this  work, 
it  is  essencial  for  clarity  that  we  briefly  explain 
the  meaning  of  the  words  "energy  crapping."  Since 
che  pure  chickness-extensional  resonant  frequencies 
are  cutoff  frequencies,  there  is  usually  a  nearby 
frequency  range  in  which  the  transverse  mode  shape 
is  evanescent.  There  is  also  a  nearby  frequency 
range  in  which  che  cransverse  mode  shape  is  trig¬ 
onometric.  Consequently,  by  che  selection  of  Che 
appropriate  chickness-extensional  overtone  (or 
fundamental)  and/or  the  appropriate  adjustment  of 
the  geometry  in  the  electroded  and  unelec  eroded 
regions,  che  cransverse  modal  behavior  can  be  made 
to  decay  with  distance  away  from  che  electrodes  in 
che  unelectroded  region.  The  resulting  vibration 
is  called  a  trapped  energy  mode,  which  radiaces  a 
concrollably  small  amount  of  energy  inco  the  adja¬ 
cent  thick  porcion  of  Che  semiconducting  wafer  and, 
hence,  results  in  che  highest  possible  Q,  albeit 
with  many  nearby  spurious  modes  with  high  Q. 
Alternatively,  che  overtone  and/or  geometry  can  be 
selected  so  that  che  mode  does  not  decay  with 
distance  away  from  the  electrode  in  the  unelec- 
croded  region  and  the  resulting  vibracion  is  not  a 
trapped  energy  mode.  In  this  case  much  more 
energy  is  radiated  into  the  adjacent  thick  porcion 
of  che  semiconducting  wafer  and  much  lower  Q's 
result.  Although  there  are  still  many  nearby 
spurious  modes,  they  are  less  troublesome  because 
the  Q's  are  lower.  All  experimental  work  on  che 
composite  resonacor  reported  to  dace  has  been  for 
this  latter  case1'5  .  On  che  other  hand  a  detailed 
analytical  treatment  of  the  composite  resonacor 
for  che  case  when  trapping  is  present  appears  in 
the  liceracure7  along  with  a  detailed  discussion 
of  when  crapping  is  and  is  not  present. 
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Using  the  existing5  constants  of  aluminum-nicride, 
we  have  found  that  the  fundamental  essentially 
chickness-extensional  mode  will  not  crap  for  an 
aluminum-nicride  film  on  a  gallium-arsenide 
diaphragm  in  che  flat  place  configuration'  .  How¬ 
ever,  the  fundamental  mode  will  crap  if  the 
gallium-arsenide  diaphragm  is  appropriately  notched 
a  small  amount  in  che  electroded  region,  as  shown 
in  Fig. 9  of  Ref. 7.  In  addition,  we  have  found  that 
the  second  essentially  thickness-extensional  mode 
will  trap  in  che  flat  plate  configuration  for  che 
same  film  and  substrate  materials  for  any  racio  of 
film-thickness  to  diaphragm- ch ickness . 

In  this  work  an  analysis  of  che  vibrations  of  a 
composite  resonacor,  which  is  driven  by  che  appli¬ 
cation  of  an  a.c.  voltage  across  scrip  eleccrodes 
on  che  major  surfaces  of  che  film,  is  performed. 

The  analysis  includes  the  pertinent  waves  in  Che 
active  region  of  the  composite  resonator,  as  well 
as  all  radiating  waves  in  che  chick  gallium- 
arsenide  plate.  All  previous  analytical  work 
expressly  ignores  radiation  inco  the  bulk  semi¬ 
conductor  except  one  treatment10,  which  unrealis¬ 
tically  ignores  the  junction  between  the  ecched 
diaphragm  and  che  bulk  semiconductor.  The  solution 
is  obcained  by  satisfying  the  differencial  equa¬ 
tions  for  the  piezoelectric  film  and  semiconductor 
as  well  as  all  boundary  conditions  on  che  major 
surfaces  of  che  film  and  semiconduc cor  exactly  and 
using  the  aepropriace  variacional  principle"  to 
satisfy  che  remaining  conditions  along  the  minor 
interfaces  approximately.  The  minor  interfaces 
separate  che  electroded  from  the  unelectroded 
regions  of  che  resonator  and  the  thin  region  of  che 
gallium-arsenide  from  the  chick  region.  Past 
experience  shows  chat  this  type  of  approximation 
yields  extremely  accurate  results  if  all  the  proper 
waves  are  included12 ,13  .  Boch  che  configuration  in 
which  the  film  ends  at  che  edges  of  che  electrodes 
and  in  which  it  continues  to  che  edges  of  the 
etched  diaphragm  have  been  considered  when  trapping 
is  not  present,  along  with  the  latter  configuration 
when  crapping  is  present.  In  each  instance  the  Q 
at  the  resonance  condition  is  calculated. 

II.  Basic  Equations 

The  differential  equations  of  motion  and  electro¬ 
statics  and  linear  piezoelectric  constitutive  equa¬ 
tions  for  che  piezoelectric  film  may  be  written  in 
the  censor  form14'15 
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TLi  - p  v  <t-°>  (l> 

„f  f  f  ^  f  f 
Tij  “  cijkXuk,X  ‘Tcij'P.k  * 

nf  f  f  f  f  ... 

V6ikA,r£ikV,  (2> 

where  the  nocation  is  conventional.  The  equations 
of  motion  and  the  linear  elastic  constitutive 
equations  for  the  gallium-arsenide  may  be  written 
in  the  censor  form15'16 

_s  s..s  _s  s  s  . 

ij,i  Uj  ’  ij  ijk2Uk,X  ' 

and  we  noce  thac  the  equations  are  too  cumbersome 
to  write  ouc  in  decail.  Ac  this  poinC  we  intro¬ 
duce  a  Cartesian  coordinate  syscem  x^,  x^  with 

the  x^-axis  normal  to  the  major  surfaces  of  the 

aluminum  nitride  film  and  along  a  cube  axis  of  Che 
gallium-arsenide.  Since  aluminum-nitride  has  small 
piezoeleccric  coupling,  for  small  wave  and  decay 
numbers  along  the  place,  which  are  the  only  ones  of 
interest  in  the  active  region,  we  need  retain  only 
x^-dependence  of  all  electrical  variables  and  we 
have 

°3 , 3  =  °  *  °3  "  ®3kXuk,i  '  £33',3  ’  (4) 


Figure  1  shows  a  schematic  diagram  of  a  cross  sec¬ 
tion  of  a  thin  aluminum-nitride  film  on  a  chin 
gallium-arsenide  layer  composite  trapped  energy 
resonator  for  the  configuration  in  which  the  film 
ends  at  che  edge  of  the  diaphragm.  In  che  uneiec- 
croded  regions  the  boundary  conditions  on  che 
major  surfaces  are 


„s  ,  ,  *..s 
i  ^  *  h  a  u 

JJ  J 


f  s 

u  .  »  u  , 
J  J 


■:  =  0  at  x,  *  0 


where  the  superscripts  f  and  s  scand  for  che  chin 
film  and  semiconductor  layer,  respectively,  and  o" 
is  che  mass  density  of  the  ground  electrode.  The 
electrical  condition  in  (3)2  is  a  consequence  of 

the  fact  chac  che  x^-  and  x2-dependence  of  all 

electrical  variables  has  been  left  out  of  account 
and  the  eleccrical  potential  in  space  is  bounded  ac 
‘<2  =  Since  the  electrodes  are  much  thinner 

chan  either  che  thin  film  or  che  layer,  we  have 
employed  approximate  thin  plate  equations1,7  for 
the  electrode  plating  in  (6)^,  and  we  have  made  use 

of  che  fact  thac  che  mechanical  stiffness  of  the 
very  chin  electrode  plating  is  negligible  for  small 
wavenumbers  along  the  plate.  Similarly,  on  the 
major  surtaces  of  che  composite  plate  in  che  elec- 
croded  region  che  boundary  conditions  are 

T  *•  ;  h  u  ,  -p  *  Ve  ac  x^  =  h  ,  (31 

•/  =  0  it  3  0  ,  [')) 


along  with  (6)  and  (7),  which  remain  unchanged. 

The  quantity  p '  is  the  mass  density  of  the  upper 
electrode. 

Figure  2  shows  a  schematic  diagram  of  a  cross- 
section  for  the  configuration  in  which  the  film 
ends  at  the  edges  of  the  electrodes.  The  boundary 
conditions  on  the  major  surfaces  of  the  composite 
plate  in  the  electroded  region,  which  still  contains 
the  film,  are  exactly  che  same  as  in  che  previous 
case.  However,  the  boundary  conditions  on  the 
major  surfaces  of  the  unelectroded  plate,  which  no 
longer  contains  the  film,  are  much  simpler  chan  in 
Che  previous  case  and  in  place  of  (5)  -  (7)  we  have 

T^  -0  at  x^  *  0  and  x^  =  -  hs  .  (10) 

For  eicher  conf i guracion  che  corresponding 
boundary  conditions  on  che  major  surfaces  in  che 
region  of  the  bulk  semiconductor  are 

T^  1  0  at  Xj  =  0  and  x^  *  -  hS  .  (11) 


Since  che  piezoelectric  coupling  is  small  in 
aluminum-nitride,  we  may  transform  the  innomo- 
geneity  from  the  boundary  conditions  in  (8),  into 


Che  differencial  equations  by  means  of  che  trans¬ 
formation  shown  in  Eqs.(3.26)  of  Ref. 7,  which  is 
considered  to  be  part  of  this  work1  s  .  Furthermore, 
since  we  have  replaced  (1)2  and  (2),  by  (4),  no 

electrical  conditions  need  be  satisfied  ac  a  minor 
interface.  Then  che  conditions  thac  should  be 
satisfied  at  each  minor  interface  are  the  conti¬ 
nuity  of  u^,  u^,  T^3  anc*  ^  .  However,  since  in 

che  approximation  technique  we  employ  che  solution 
is  written  as  a  sum  of  eigensolutions ,  each  of 
which  satisfies  the  homogeneous  differential  equa¬ 
tions  and  boundary  conditions  on  che  major  surfaces 
for  one  of  the  three  regions,  i.e.,  che  electroded 
composite,  the  unelectroded  composite  and  the  bulk 
semiconductor,  exaccly,  we  cannot  sacisfy  the 
continuity  conditions  across  the  interfaces  beeveen 
the  regions  exaccly.  Nevertheless,  the  remaining 
continuity  conditions  across  the  interfaces  may  be 
satisfied  approximately  along  with  the  inhomogen¬ 
eous  driving  cerm  resulting  from  che  aforementioned 
transformation  by  satisfying  che  proper  form  of  che 
appropriate  variational  principle’"1  of  linear 
piezoelectricity,  in  which  all  thac  remains  is  an 
integral  over  che  inhomogeneous  forcing  cerm  plus 
integrals  over  che  interfaces  because  all  other 
terms  in  that  form  of  the  variational  principle 
vanish  on  account  of  che  equations  and  conditions 
satisfied  by  the  solution  functions  emploved"  =  . 

The  form  of  che  variacional  principle  of  interest 
here  is  given  in  Eq.(6.44)  of  Ref.  11.  Since  the 
solution  functions  sacisfy  the  aforementioned 
equations  and  conditions,  all  chac  remains  of 


Eq. (6.44)  of  Ref . 11  is 


dx^dx 
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Similarly,  for  the  completely  electroded  composite 
plate  driven  by  a  voltage  it  has  been  shown  that 
the  pure  thickness-extensional  resonant  frequencies 
are  given  by 

(cX3/cf)1/2T°(l+P0/G0)  ,  (IS) 


_  S  \  /  c  s  ,»S,  ,  S  *S,  _S  ,  - Z,s  s 

T13H6U3  +  6u^)  +  (ux  -  u1)(£T11  +  5T11) 

+  (u3-G3)(5T13  +  JX13)]x1=ddx3*0'  <12) 

for  the  case  when  trapping  is  noc  present  and  the 
film  extends  over  che  entire  diaphragm.  Mien  the 
film  ends  at  the  edges  of  Che  electrodes  the  fifth 
incegral  in  (12)  is  not  present  and  the  second 
integral  is  replaced  by 


J(?i>3f  +  xn5u'if)l  .  dV  (13) 
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k-»  (e‘3)*/e*3e33,  R'  =  :'h7;h\  (1?) 

The  difference  between  the  resonant  frequencies  of 
pure  thickness-extensional  vibrations  of  che 
composite  plate  vich  and  without  a  driving  elec¬ 
trode.  i.e.,  x  -  -u  for  the  same  -J,.  is  very 
’  e  e  f 

important  when  trapping  is  present. 


When  trapping  is  present  wich  a  nocch  under  che 
electroded  region,  as  shown  in  Fig. 3,  the  fourth 

incegral  is  taken  from  -hS  to  0  (instead  of  -hS 
to  0)  and  there  is  an  additional  terra 


fTn5uI  +Tn5u3)  ,  dx3 

V1 


in  Eq. (12). 


IV.  Straight-Crested  Dispersion  Relations 

In  this  section  we  obtain  the  straight-crested 
dispersion  relations  for  the  composite  plate,  che 
thin  diaphragm  without  the  film  and  che  bulk  semi¬ 
conductor  because  these  determine  the  solution 
functions  that  are  used  in  che  variational  condi¬ 
tion  (12).  Since  these  eigensolutions  in  each 
region  sacisfy  the  homogeneous  differential  equa¬ 
tions  and  boundary  conditions  on  che  major  surfaces, 
they  are  coo  cumbersome  to  present  here'5. 


III.  Pure  Thickness-Extensional  Vibrations 

Since  in  all  cases  the  composite  resonator  will  be 
operated  at  a  frequency  in  the  vicinity  of  one  of 
the  pure  thickness-extensional  resonances  of  the 
composite  place,  the  pure  thickness-extensional 
resonant  frequencies  of  che  composite  are  of 
particular  importance  in  this  work.  It  is  shown 
in  Sec. Ill  of  Ref. 7  that  for  the  composite  place 
without  a  driving  electrode  but  with  a  ground 
electrode,  che  pure  thickness-extensional  resonant 
frequencies  are  given  by 

"e*  <4/‘’V/2H°f(l  -  R"/G°)  ,  (15) 

o 

where  is  the  appropriate  root  of 


The  solution  for  decaying  waves  in  the  unelectroded 
composite  plate  is  presented  in  Eqs.(4.14)-  (4.24) 
of  Ref. 7,  which  is  considered  to  be  part  of  this 
work.  As  noted  in  Ref. 7,  the  solutions  for  travel¬ 
ing  waves  in  the  unelectroded  and  electroded 
composite  place  can  readily  be  obtained  from  the 
solution  presented.  The  solutions  for  the  simpler 
cases  of  the  diaphragm  without  the  film  and  the 
bulk  semiconductor  are  not  as  cumbersome,  but  are 
also  too  lengthy  to  present  here'5.  These  types 
of  solutions  are  presented  in  a  number  of 
places1  s  ,19  . 

The  pertinent  dispersion  curves  for  the  aluminum- 
nitride  film  on  the  gallium-arsenide  layer  compo¬ 
site  plate  are  shown  in  Fig. 4.  We  do  not  distin¬ 
guish  beeween  the  electroded  film  wich  shorted 
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electrodes  and  the  unelectroded  film  In  the  figure. 
This  is  done  in  greac  decail  in  Ref. 7.  Further¬ 
more,  when  required, Che  difference  between  the  two 

g  g 

can  be  calculated  from  (<u  -  tu  )  given  in  the  last 
section.  The  dispersion  curves  are  for  a  film 
thickness  of  7  microns  and  a  diaphragm  thickness 
of  14  microns.  The  pertinent  dispersion  curves 
for  a  14  micron  chick  gallium-arsenide  diaphragm 
are  shown  in  Fig. 5.  The  lowest  17  real  dispersion 
curves  for  the  bulk  gallium-arsenide  plate  are 
shown  in  dimensionless  form  in  Fig. 6.  For  the 
fundamental  essentially  thickness-extensional 
resonance  for  the  composite  dimensions  mentioned 
(around  132  MHz),  this  number  of  dispersion  curves 
is  for  a  gallium-arsenide  wafer  thickness  of  about 
5  rails.  In  chis  work  we  perform  calculations  for 
wafer  thicknesses  up  Co  8  mils  for  which  there  are 
30  real  dispersion  curves  for  a  frequency  of 
132  MHz.  However,  we  do  not  bother  to  show  the 
figure  for  more  chan  17  in  chis  work. 

V.  Forced  Vibrations  of  Composite  Resonator 

In  chis  section  we  determine  the  essentially 
thickness-excensional  vibrations  driven  by  Che 
application  of  a  sceady-scate  driving  voltage  to 
the  scrip  electrodes  on  Che  surfaces  of  the  piezo¬ 
electric  film  of  the  composite  resonator  shown  in 
any  of  Figs.  1-3.  Since  we  include  radiation  inco 
Che  bulk  semiconductor,  we  can  use  the  solution  to 
calculate  Che  Q  at  resonance  of  the  mode  resulting 
from  radiacion  inco  the  bulk  semiconductor. 


In  accordance  with  the  earlier  discussion  we  cake 
:ne  approximate  solution  in  the  torn 
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where  h,  n  and  n  are  Che  number  of  branches  of  che 
dispersion  curves  included  in  che  eleccroded 
composite,  unelectroded  composite  or  diaphragm  and 
bulk  semiconductor,  respectively,  which  are  given 

in  Figs.4-6.  The  eigensolucion  functions  u'^, 

-s(3>  i(-j)  sta)  .s  .  L  .  3 

u  ,  u  ,  u  ,  u  denote  che  solution  tunc- 
a  a  a  ’  a 

cions  chat  satisfy  the  differential  equations  and 
die  boundary  conditions  on  the  major  surfaces  tor 
each  of  Che  respective  regions  and  yield  che  dis¬ 
persion  curves  shown  in  Figs. 4-  6  and  which  are 
coo  cumbersome  to  present  here' 3 .  Since  the  solu- 

—  (  a  > 

cion  tunc  cions  in  (30)  are  fixed,  only  the  K'"  , 

and  ^  are  varied  when  (20)  is  substituted 
into  ill).  Accordingly,  substituting  from  (20' 


into  (12),  employing  (2)^  and  (3)j  and  performing 
the  integrations,  we  obtain15  an  equation  of  che 
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where  the  expressions  for  c'  .  a,  .  a  .a  b 

’  35’  o e’  vC’  3o 

and  b_^  are  coo  lengthy  to  present  here.  Since 
the  variations  in  (21)  are  arbitrary,  we  obtain' f 

—  /  ft  \ 

the  inhomogeneous  linear  equations  for  the  ^ 
and  in  the  form 
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VI.  Quality  Factor  Resulting  from  Radiation 

In  chis  section  we  calculate  che  quality  factor  due 
to  radiation  inco  the  bulk  semiconductor  for  each 
of  the  three  cases  from  the  analysis  presented  in 
che  previous  section.  Although  the  solution  may 
be  obcained  from  the  analysis  in  che  previous  sec¬ 
tion  at  any  driving  frequency  for  which  all  che 
pertinent  dispersion  curves  are  available,  in  this 
work  we  are  interested  in  the  solution  only  ac  the 
fundamental  essentially  thickness-excensional 
resonance  for  each  of  the  chree  configurations. 

It  is  clear  from  pasc  experience' ~  3  chat  under 

these  circumstances  we  need  consider  only  the 
essentially  thickness-excensional  branch  in  the 
composite  region  (either  eleccroded  or  noc)  of  the 
resonator,  which  is  che  curve  labeled  1  in  Fig. 4. 
However,  in  che  chin  region  of  the  semiconduc cor 
withouc  che  tilm  3nd  che  bulk  semiconductor  all 
pertinenc  dispersion  curves  shown  in  che  respective 
Figs. 5  and  0  muse  be  included  co  obtain  accuracv. 
The  use  of  only  one  branch  for  the  composite  region 
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means  that  we  always  have  n«  1  and  n»  1  for  the 
configuration  shown  in  Figs.l  and  3  but  not  for  the 
configuration  shown  in  Fig. 2,  for  which  all  curves 
shown  in  Fig. 5  must  be  included. 

As  usual,  the  quality  factor  Q  is  defined  by 

Q-(K  +  U)/ECU,  (23) 


4.  VJ 

K  ’  T  J  dC  2  J  dXl  I  2  °VadX3> 


u  =  t  I  dc  :  Re[  I  dxi  J  2  (cnui,iul,i 

o  -d  -hS 

+  C33U3,3U3,3  +  2C13U1,1U*3,3+C44(U1J3 

+  u3,l)(u1.3+U3,ir)dX3]  > 

T  ’  o 

EtW  *  -  2  f  dt  -g  Re  f  (-T.  u~ )  ,  dx  (24) 

J  2  J  la  a  x  *d  3 


in  which  T  is  the  period  of  Che  vibration  and  we 
have  taken  the  liberty  of  writing  the  integrals  in 
(24)^  ,,  over  discontinuous  functions"5  to  achieve 

brevity.  For  a  given  geometry  and  mode,  resonance 
is  determined  by  obtaining  Q  over  a  (narrow)  range 
of  frequencies  and  finding  the  frequency  for  which 
Q  is  a  maximum. 


In  performing  che  calculations  we  have  found  that 
it  is  imperative  that  all  radiating  plate  waves  in 
the  thick  region  of  che  gallium-arsenide  be 
included  in  order  to  achieve  accuracy.  Since  at  a 
given  frequency  the  number  of  radiating  waves  in  a 
plate  goes  up  significantly  with  thickness,  we  have 
considered  gallium-arsenide  wafers  no  thicker  than 
8  mils  at  a  frequency  around  132  MHZ,  for  which 
there  are  30  radiating  place  waves.  Results  are 
presented  for  wafer  thicknesses  ranging  from  4  mils 
to  8  mils  because  this  is  considered  to  be  within 
che  practical  range.  All  che  results  presented  are 
for  a  film  thickness  of  7  microns  and  a  diaphragn 
thickness  of  14  microns  and  che  laceral  dimensions 
of  each  configuration  were  adjusted  slightly  to 
maintain  the  same  resonant  frequency  for  computa¬ 
tional  convenience.  Since  the  calculated  Q  is  a 
very  rapidly  varying  function  of  the  wafer 
thickness .calculations  had  to  be  performed  for 
very  small  increments  in  thickness,  i.e.,  1  micron, 
in  order  to  get  all  the  peaks  and  valleys  in  the 
interval . 


in  Q  with  wafer  thickness  similar  to  Fig. 7,  but 
in  this  case  the  Q's  are  considerably  higher.  It 
can  be  seen  from  the  figure  that  the  highest  Q 
obtained  is  about  15,000  and  the  lowest  Q  is  about 
200.  Calculations  were  performed  when  trapping  was 
induced  in  the  fundamental  mode  by  notching  che 
diaphragm  under  che  electrode  as  shown  in  Fig. 3. 

As  expected,  the  Q  due  to  radiacion  can  be  made  as 
high  as  we  wish  simply  by  extending  che  laceral 
dimensions  of  the  film  and  diaphragm.  Figure  9 
gives  Q  as  a  function  of  (d  -  l) ,  i.e.,  the  distance 
from  the  edge  of  the  electrode  to  the  edge  of  the 
diaphragm.  It  can  be  seen  from  the  figure  that  Q 
increases  very  rapidly  with  (d  -  L)  and  at  d  -  21, 

Q  *  70,000  for  this  geometry.  For  these  calcula¬ 
tions  the  wafer  thickness  was  6  mils.  Even  when 
trapping  is  present  che  Q  is  a  varying  function  of 
uafer  thickness,  but  the  range  of  the  variation  is 
less  chan  1/10  of  chat  when  trapping  is  not  present, 
which  is  not  of  interest  since  the  Q  with  crapping 
is  so  much  higher  chan  che  Q  without  trapping.  In 
interpreting  the  foregoing  information  it  should  be 
remembered  that  the  high  Q's  calculated  should  be 
higher  than  che  accual  Q's  because  the  material  Q 
and  che  Q  due  to  radiation  into  the  air  are  not 
included.  The  results  discussed  above  mean  that  in 
order  to  obtain  reasonably  high  Q  when  trapping  is 
not  present, for  a  given  wafer  thickness  the  thick¬ 
nesses  of  the  film  and  diaphragm  must  be  very 
precisely  selected. 

Since  the  variational  condition  in  (12)  does  not 
expressly  match  che  mode  shape  at  the  interfaces 
unless  the  solution  is  exact*5,  the  extent  to  which 
the  calculated  mode  shape  matches  at  the  inter¬ 
faces  gives  an  indication  of  the  accuracy  of  the 
approximate  solution  obtained.  Typical  plots  of 
the  u^-displacenent  field,  which  is  che  large  one 

for  the  essentially  thickness-extensional  modes 
considered  here,  for  the  configurations  shown  in 
Figs.l  and  2  are  shown  in  Figs. 10  and  11.  respect¬ 
ively.  It  can  be  seen  from  the  figures  chat  che 
Uj-displaceroent  field  macches  quite  well  at  the 
interfaces.  Consequently,  we  can  conclude  thac  the 
approximate  solution  obcained  is  quite  accurate. 
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Figure  l  Cross-Section  of  a  Composite  Resonacor 
Consisting  of  a  Thin  Piezoelectric  Film 
on  a  Gallium-Arsenide  Layer  with  t-.e  Film 
Covering  the  Encire  Layer 


Figure  2  Cross-Section  of  a  Composite  Resonator 
wich  the  Piezoelectric  Film  Ending  at 
the  Edge  of  the  Electrodes 
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Figure  3  Cross-Section  of  a  Composite  Resonator 
with  the  Layer  Notched  under  the  Elec¬ 
trodes  to  Cause  Trapping  in  the  Funda¬ 
mental  Mode  when  it  does  not  Trap  in  the 
Flat  Plate  Configuration 


Figure  L  Dispersion  Curves  for  the  Composite 

Region  Consisting  of  an  Aluminum-Nitride 
Film  on  a  Gallium-Arsenide  Layer  wich  5 
she  Propagation  Wavenumber 


Figure  5  Dispersion  Curves  for  the  Gallium- 

Arsenide  Layer  without  the  Thin  Film 


Figure  6  Dispersion  Curves  for  the  Callium- 

Arsenide  Wafer  with  Z  the  Dimensionless 
Frequency  Normalized  with  Respect  to  the 
First  Thickness-Shear  Frequency. 


Figure  7  Quality  Factor  Versus  Wafer  Thickness 
when  Trapping  is  not  Present  for  the 
Composite  Resonator  Configuration 
Shown  in  Figure  1 


Figure  8  Quality  Factor  Versus  Wafer  Thickness 
when  Trapping  is  not  Present  for  the 
Composite  Resonator  Configuration 
Shown  in  Figure  2 
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Quality  factor  of  the  piezoelectric  thin  film  on  semiconductor  composite 
resonator  resulting  from  radiation  into  the  semiconductor  wafer 

D.  V.  Shick,  D  S.  Stevens,'"  and  H.  F.  Tiersten 
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( Received  3  March  ll>86;  accepted  for  publication  7  May  14S6) 

The  compoMtc  resonator  consists  of  a  uniform  thin  layer  etched  in  a  small  well-defined  legion 
of  a  semiconductor  wafer  to  form  a  diaphragm,  upon  which  is  deposited  a  thin  piezoelectric 
ft  Im  along  with  the  electrodes  to  form  a  resonant  region  directly  on  I  he  wafer  Although  the 
composilo  resonator,  which  operates  m  an  essentially  thickness-extensional  mode,  can  he 
c*  oisi  rncicu  :•  >  employ  energy  trapping,  almost  all  existing  experimental  work  in  the  literal  me 
is  for  I  lie  case  w  hern  napping  is  not  present  All  prev  ions  analytical  work  expresslv  ignores 
rad'.itK  n  mio  the  hulk  semiconductor  except  one  treatment,  which  unrealixticalh  ignores  the 
pm. lion  between  I  he  etched  diaphragm  and  the  bulk  semiconductor  In  this  work  an  analysis 
of  the  compo  ilc  resonator  driven  into  essentially  thickness-extensional  vibrations  by  the 
application  of  a  voltage  to  strip  electrodes  is  performed.  The  analysis  includes  all  radiating 
plate  waves  in  the  thick  portion  oft  he  semiconductor.  The  solution  consists  of  a  sum  of  terms 
satisfying  all  differential  equations  and  boundary  conditions  on  major  surfaces  exactly  and  uses 
the  appropriate  variational  principle  of  linear  piezoelectricity  to  satisfy  the  remaining 
conditions  approximately,  l  or  the  case  of  the  aluminum  nitride  film  on  gallium  arsenide  the  Q 
is  calculated  for  both  the  configuration  in  which  the  film  ends  at  the  edges  of  the  electrodes 
and  in  which  it  continues  to  the  edges  of  the  etched  diaphragm  when  trapping  is  not  present 
and  for  the  latter  configuration  when  trapping  is  present.  The  calculations  show  that  when 
trapping  is  not  present  the  Q  is  a  very  rapidly  vary  ing  function  of  the  ratio  of  the  composite 
resonator  thickness  to  the  wafer  thickness  and  that  the  range  of  variation  is  very  large,  i.e., 
between  one  and  two  orders  of  magnitude.  The  calculations  also  reveal  that  when  trapping  is 
present  the  Q  is  alw  ay  s  much  larger  and  its  range  of  variation  much  smaller  than  when 
trapping  is  not  present. 


I.  INTRODUCTION 

The  composite  resonator  consists  of  a  uniform  thin  layer 
etched  in  a  small  well-defined  region  of  a  semiconducting 
wafer  to  form  a  diaphragm,  upon  which  is  deposited  a  thin 
piezoelectric  film  along  with  the  electrodes  to  form  a  reso¬ 
nant  region  directly  on  the  wafer.  In  this  work  an  analysis  of 
the  vibrations  of  a  composite  resonator  is  presented  includ¬ 
ing  the  very  important  transverse  behavior  of  the  essentially 
thickness-extensional  modes  and  the  attendant  radiation 
into  the  bulk  semiconductor  both  for  cases  in  which  energy 
trapping  is  and  is  not  present.  The  calculations  are  per¬ 
formed  for  the  particular  ease  of  an  aluminum-nitride  film 
on  gallium  arsenide. 

before  proceeding  w  ith  a  discussion  of  this  work,  it  is 
essential  for  clarity  that  we  briefly  explain  the  meaning  of  the 
words  "energy  trapping."  Since  the  pure  thickness-exten¬ 
sional  resonant  frequencies  are  cutoff  frequencies,  there  is 
usually  a  nearby  frequency  range  in  which  the  transverse 
mode  shape  is  evanescent  There  is  also  a  nearby  frequency 
range  in  which  the  transverse  mode  shape  is  trigonometric. 
Consequently .  by  the  selection  of  the  appropriate  thickness- 
extensional  overtone  for  fundamental )  and/or  tile  appro¬ 
priate  adjustment  <  >f  tie  geometry  in  the  elect  roiled  and  une- 
leclroded  lemons,  the  iiansverse  modal  behavior  can  be 
mad-  lodeeav  w  it  li  dl- lance  aw  ay  from  I  he  elect  lodes  ill  I  he 
uncle  ■to  .led  legion  I  he  icstiltinc  vibialtoti  is  called  a 
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trapped  energy  mode,  which  radiates  a  controllable  small 
amount  of  energy  into  the  adjacent  thick  portion  of  the  semi¬ 
conducting  wafer  and.  hence,  results  in  the  highest  possible 
Q,  albeit  with  many  nearby  spurious  modes  with  high  Q 
Alternatively .  the  overtone  and/or  geometry  can  be  selected 
so  that  the  mode  does  not  decay  with  distance  away  from  the 
electrode  in  the  unelect roded  region  and  the  resulting  vibra¬ 
tion  is  not  a  trapped  energy  mode.  In  this  case  much  more- 
energy  is  radiated  into  the  adjacent  thick  portion  oft  he  semi¬ 
conducting  wafer  and  much  lower  Q  \  result.  Although 
there  are  still  many  nearby  spurious  modes,  they  are  less 
troublesome  because  the  Q\  are  lower.  Almost  all  experi¬ 
mental  work  on  the  composite  resonantor  reported  to  date 
has  been  for  this  latter  case  1  *’  On  the  other  hand,  a  detailed 
analytical  treatment  of  the  composite  resonator  for  the  case 
when  trapping  is  present  appears  in  the  literature  along  w  it  li 
a  detailed  discussion  of  when  trapping  is  and  is  not  present 
Using  the  existing'  constants  of  aluminum  nitride,  we 
have  found  that  the  fundamental  essentially  thickness-ex- 
tensional  mode  will  not  trap  for  an  aluminum-nitride  film  on 
a  gallium  arsenide  diaphragm  m  the  flat-plate  configura¬ 
tion  '  However,  the  fuiulamcnl.il  mode  will  trap  it  The  gal¬ 
lium  arsenide  diaphragm  o  uppiopriutcly  notched  sm.J1 
amount  m  the  elect  roded  u  cion,  as  show  n  in  I  le  '<  ■  t  Rl  I 
In  addition,  we  have  found  that  'he  second  cssenll..  v  i  ■ 1 1.  k 
ness-cxtcnsional  mode  will  itap  m  the  flat -plate  e  :i!..-u:  i 

lion  for  llu  same  film  and  Mibo  :  ,,h  uiutcri  i!s  h  a  ..."  : 1 
film  thickness  to  diaplll  agm  t  l.n  Kuess 
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resonator,  which  is  driven  by  the  application  of  an  ac  voltage  onator  and  the  thin  region  of  the  gallium- arsenide  from  the 

across  strip  electrodes  on  the  major  surfaces  of  the  film,  is  thick  region.  Past  experience  shows  that  this  type  ofapproxt- 

performed.  The  analysis  includes  the  pertinent  waves  in  the  mat  ion  yields  extremely  accurate  results  if  all  tlu  proper 

active  region  of  the  composite  resonator,  as  well  as  all  radiat-  waves  are  included. 1 '  1 '  Both  the  configuration  m  w  Inch  the 

ing  waves  in  the  thick  gallium-arsenide  plate.  All  previous  film  ends  at  the  edges  of  the  electrodes  and  in  which  u  con- 

analytical  work  expressly  ignores  radiation  into  the  bulk  se-  tinues  to  the  edges  of  the  etched  diaphragm  have  been  con- 

iniconductor  except  one  treatment,1"  which  unrealistically  sidered  when  trapping  is  not  present,  along  with  the  latter 

ignores  the  junction  between  the  etched  diaphragm  and  the  configuration  when  trapping  is  present.  In  each  instance  the 

bulk  semiconductor.  The  solution  is  obtained  by  satisfying  quality  factor  ( Q )  at  the  resonance  condition  is  calculated, 

t he  differential  equations  for  the  piezoelectric  film  and  semi-  The  calculations  show  that  w  hen  trapping  is  noi  present  t he 

conductor  as  well  as  all  boundary  conditions  on  the  major  Q  is  a  very  rapidly  vary  ing  function  of  the  i  .1110  of  the  corn- 

surfaces  of  the  film  and  semiconductor  exactly  and  using  the  posite  resonator  thickness,  t.c  .  the  piezoelectric  til  in  plus  the 

appropriate  variational  principle."  in  which  all  conditions  semiconductor  layer,  to  the  wafer  thickness  and  that  the 

are  natural  conditions,  to  satisfy  the  remaining  conditions  range  of  variation  of(Jis  very  large,  i  c  .  between  one  and  two 

along  the  minor  interfaces  and  the  inhomogeneous  driving  orders  of  magnitude.  The  calculations  also  u  -cal  that  when 
term  that  appears  in  a  differential  equation  for  the  film  in  the  trapping  is  present  the  Q  is  alw  ay  s  much  !  irgei  and  its  .  .nice 

clectroded  region  approximately  The  minor  interfaces  sepa-  of  variation  with  the  thickness  rain  \  much  smaller  than  w  hen 

rate  the  eleetroded  from  the  unelectroded  regions  of  the  res-  trapping  o  not  present 

- 1 

II.  BASIC  EQUATIONS 

The  differential  equations  of  motion  and  electrostatics  and  linear  pie/oelectiic  constitutive  equations  fm  aluminum 
nitride  with  .v,  the  hexagonal  axis  may  he  written  in  the  form 1 1 
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w  heie  >,  /  .  arid  !>  denote  the  components  of 'lie  mechanical  displacement,  the  stress  lensoi .  and  the  i  lectnc  displacement, 

respective! ,.  and  <  <•  .  atul  r  denote  the  elastic,  piezoelectric,  and  didcviiic  constants,  t .  spi  vttvdy .  y  and  /•  denote  ills 

elect rtc  potential  and  muss  density .  1  espccliv  cly .  and  we  have  employ  ed  the  usual  compressed  matrix  notation  \'  e  ctupf  t> 
(  a;  test  an  i  .tisor  total  ion.  the  summation  eon  ,  eniion  for  i  cpi  ated  tensor  indices,  the  dot  notation  for  partial  dill  etc  to  i.ition 
'■vi1  h  respes  i  to  lime,  and  t  he-  convent  n  -n  that  a  comma  followed  by  an  index  denotes  pari  ml  differentiation  with  respect  in  a 
pac  'em  i due tc  At  this  point  we  note  that  w e  have  int reduced  a  (.  u tesian  coot d male  sv stem  v ,,  v  with  the  v ,  axis  normal 
to  the  m.iioi  surfaces  ol  the  aluminum  nitride  film.  Similarly,  the  displacement  equations  of  motion  and  linear  clastic 
sonsfi’utne  equations  for  the  nonpiezodectric  gallium  arsenide  semiconductor  with  v .  along  a  cube  avis  lake  the  form" 
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A  a-  -w  note  that  in  the  essentially  thiekness-extension- 
d  i-i' d  s  tn  the  region  containing  the  film,  which  are  the 
'  ".tsiest  rn  tho  work,  the  wave,  or  decay  mtmhets 
ai-  -  ..  ;  plat  .-  are  much  smaller  film  the  ihtekness  wave- 

i  i  it  or-  X-  i  i.  onse'mence  ol  these  small  wave  and  dec  p 
•  m-tb  •  w  the  |  -l;*t%-  amt  the  small  pH,  oeleclr  ic  coupling 

'■  cm  mum  nitride,  we  cm  ignore  the  v  and  \  depen 
■  nee  nt  ah  i  lee' rival  v.uial'1-  s  ,is  m  Ret  '.  and  in  place  ol 
I  ,qs  i  2  i  and  <  4 } .  we  have 


c ,  ,u  .  <  .  tt.  .rn 

I)  e.  u  <  i  1  'i 

I  igui e  !  show s  a  schematic  -  I lag i  an  >t  ,i  cross  sc.  lion  ol 
a  thin  aluminum  nitride  film  on  ■  thin  gallium  aiscntde  l.o 
er  composite  tr  j  pcd-cncrgy  tvsonatoi  toi  ihc  contigmatn  n 
m  w  Inc  1 1  ti  c  til m  ends  at  the  nice  ol  the  diaphragm  In  the 
uncli  c troiled  regions  the  boundat  y  conditions  on  the  m  oot 
surfaces  am 
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MG  !.  Cross  section  «»t  a  composite  resonator  consisting  of  a  thin  piezoe¬ 
lectric  film  on  a  gallium  arsenide  layer  vulh  the  film  covering  the  entire 
lav  or 


F  ICi.  2  Cross  section  of  a  composite  resonator  with  the  piezoelectric  film 
ending  at  the  edge  of  the  electrodes. 


0  at  .v  < 


D  \  --  0  at  -V ,  -  0.  (  10) 


0  at  ,v , 


where  the  superscripts  /  and  ?  stand  for  the  thin-film  and 
semiconductor  laser,  respectively,  and/)"  is  the  mass  density 
of  the  ground  electrode.  The  electrical  condition  in  (9),  is  a 
consequence  of  the  fact  that  the  r,  and  v  ,  dependence  of  all 
electrical  variables  has  been  left  out  of  account  and  the  elec¬ 
tric  lield  in  space  vanishes  at  .v,  —  x  .  Since  the  electrodes 
are  much  thinner  than  either  the  thin  film  or  the  layer,  we 
have  employed  approximate  thill-plate  equations'  for  the 
electrode  plating  ttt  (  10) and  we  have  made  use  of  the  fact 
that  the  mechanical  stillness  of  the  very  thin  electrode  plat¬ 
ing  is  negligible  for  small  wavenumbers  along  the  plate.  Sim¬ 
ilarly.  on  the  ma|or  surfaces  of  the  composite  plate  in  the 
elect roded  ration  the  h  umdarv  conditions  are 


along  with  i  It)  i  and  (II).  which  remain  unchanged  The 
quantity  ,>  is  i he  mass  density  of  the  upper  electrode. 

f  igure  1  'Iv'iis  a  schematic  diagram  of  a  cross  section 
fa  the  eontigiu.ition  m  which  the  film  etuis  at  the  edges  of 
'he  electrodes  I  he  boundary  conditions  on  the  major  sur¬ 
faces  of  the  composite  plate  in  the  elect  roded  region,  which 
still  contains  the  lilm.  are  exactly  the  same  as  in  the  previous 
case  However,  the  boundary  conditions  on  the  major  sur¬ 
taxes  of  the  unelectroded  plate,  which  no  longer  contains  the 
film,  are  much  simpler  than  in  the  previous  case  and  in  place 
of  (9)-  (  1  1  )  we  have 


7 =  0  at  -v,  =  0  and  xx  =  ~h'.  (14) 

For  either  configuration  the  corresponding  boundary  condi¬ 
tions  on  the  major  surfaces  in  the  region  of  the  bulk  semicon¬ 
ductor  are 

/ T,  —  0  at  .v,  —  0  and  .v,  -  -  h  (1?) 

Since  we  have  replaced  (2)  and  (4)  by  (7)  and  (8).  respec¬ 
tively.  no  electrical  conditions  need  be  satisfied  at  a  minor 
interface.  Then  the  conditions  that  should  be  satisfied  at 
each  minor  interface  are  the  continuity  of  a,,  a,.  I  j ,.  and  7j  , 
and  at  each  free  edge  are  7j ,  =  7j ,  =  0.  However,  since  in 
the  approximation  technique  we  employ  the  solution  is  writ¬ 
ten  as  a  sum  of  eigensolutions.  each  of  which  satisfies  the 
homogeneous  differential  equations  and  boundary  condi¬ 
tions  on  the  major  surfaces  for  one  of  the  three  regions,  i.e.. 
the  eleetroded  composite,  the  unelectroded  composite,  and 
the  bulk  semiconductor,  exactly,  we  cannot  satisfy  the  con¬ 
tinuity  conditions  across  the  interfaces  between  the  regions 
exactly.  Nevertheless,  the  remaining  continuity  conditions 
across  the  interfaces  may  be  satisfied  approximately  by  satis- 
fving  the  proper  form  of  the  appropriate  variational  princi¬ 
ple"  of  linear  piezoelectricity,  in  which  all  conditions  are 
natural  conditions  and  all  that  remains  is  an  integral  over  an 
inhomogeneous  forcing  term  that  arises  in  the  solution  plus 
integrals  over  the  interfaces,  because  all  other  terms  in  thai 
form  of  the  variational  principle  vanish  on  account  of  the 
equations  and  conditions  satisfied  by  the  solution  functions 
employed.  The  form  of  the  variational  principle  of  interest 
here  is  given  in  Eq.  (6.44)  of  Ref.  1 1,  which  we  reproduce 
here  for  the  configuration  shown  in  Fig.  3  in  the  form"' 
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FIG.  3.  Diagram  of  a  bounded  region  containing  an  internal  surface  of  dis¬ 
continuity. 
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where  S""  denotes  the  surface  of  discontinuity  separating 
region  1  from  region  2,  and  S'”'  and  5jr'"1  denote  the  por¬ 
tions  of  the  with  surface  on  which  natural-  and  constraint- 
type  conditions,1"'  respectively,  are  prescribed.  We  note  that 
S'""  and  S\""  refer  to  different  portions  of  the  surface  for 
different  terms  in  the  boundary  sums  depending  on  each  ac¬ 
tual  condition  at  a  point.  In  ( 16)  n[m>  denotes  the  outwardly 
directed  unit  normal  to  the  with  surface,  rt[d'  denotes  the  unit 
normal  to  the  surface  of  discontinuity  directed  from  region  1 
to  region  2,  t  a  u\'"\  and  <p‘""  denote  the  prescribed 
tractions,  charge,  mechanical  displacement,  and  electric  po¬ 
tential,  respectively,  for  the  W7th  surface,  and  the  meaning  of 
the  remaining  quantities  in  (16)  is  obvious  from  earlier  dis¬ 
cussion.  The  variational  condition  in  ( 16)  for  the  simplest 
configuration,  which  is  shown  in  Fig.  3,  makes  clear  how  to 
apply  the  variational  condition  to  more  complicated  config¬ 
urations  such  as  shown  in  Figs.  1  and  2.  Consequently,  it  is 
not  worthwhile  writing  the  general  variational  equivalent  of 
(16)  for  each  configuration  because  it  becomes  too  cumber¬ 
some  and  the  actual  useful  condition  can  be  obtained  from 
(16)  for  any  configuration.  Furthermore,  in  any  given  appli¬ 
cation  what  is  taken  as  a  surface  of  discontinuity  for  applica¬ 
tion  of  the  variational  principle  (16)  is  not  clear  from  the 
figure  itself,  but  depends  on  what  conditions  the  approxi¬ 
mating  functions  satisfy. 

III.  PURE  THICKNESS-EXTENSION AL  VIBRATIONS 

Since  in  all  the  cases  the  composite  resonator  will  be 
operated  at  a  frequency  in  the  vicinity  of  one  of  the  pure 
thickness-extensional  resonances  of  the  composite  plate,  the 
pure  thickness-extensional  resonant  frequencies  of  the  com¬ 
posite  are  of  particular  importance  in  this  w  ork.  It  is  shown 
in  Sec.  Ill  of  Ref.  7  that  for  the  composite  plate  without  a 
driving  electrode  but  with  a  ground  electrode,  the  pure 
thickness-extensional  resonant  frequencies  arc  given  by 

o>.  W\,/pV  V,’<1  R"/G"),  (17) 

where  >/“  is  the  appropriate  root  of 

tan  i/./t  '  ‘  c'fi  tan  fia  tph  '  0,  (18) 

and 


R  " 


r"i> 
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cos  ?/',7t  ' 


~£L‘UL —  .  (19) 

cos-  n<r>f]h  ' 


w  it  h 


<  ,  .  '  *  |  («\,  )/<.',  ]  .  €•' 
ft  i (  p'/c\,  ft')'  ,  a  h  /h  1 . 


(20) 


Similarly,  for  the  completely  electroded  composite 'plate 
driven  by  a  voltage  it  has  also  been  shown  in  Sec.  Ill  of  Ref.  7 
that  the  pure  thickness-extensional  resonant  frequencies  are 
given  by 

5,  -  1  +P°/G °).  (21) 

where 

p<>  _  ft  ~ 

W}h')2 

X  ( - — -  —2  +  cr/u  tan  pa  tf'.h  '  tan  rf)h 

V  cos  t fjh  1  ) 

—  R  "  —  R  '(  1  —  crp  tan  rf]h  '  tan  pa  r]'}h  f),  (22) 

and 

ft 3  =  (eU  )7c{,e{,,  R  '  =  p'h  '/p’h  >.  (23) 

The  difference  between  the  resonant  frequencies  of  pure 
thickness-extensional  vibrations  of  the  composite  plate  with 
and  without  a  driving  electrode,  i.e.,  at,.  —  to,.  for  the  same 
if,,  is  very  important  when  trapping  is  present. 

For  later  use  in  this  work  we  now  note  from  Eqs.  (3.13) 
and  (3.14)  of  Ref.  7  that  when  there  is  no  driving  electrode 
the  solution  for  pure  thickness  vibrations  takes  the  form 

u{  =  ( A  'cos  i],  x,  +  B  'sin  r]f  xs)e"",  (24) 

u\  =  (A  '  cos  7;,x,  +  B  '  sin  rjsxx)e“",  (25) 

where  from  Eq.  (3.8)  of  Ref.  7  we  have 

< f  1  =  (e,', /<?,', )  [«,'  -  u{ (0,1 )  ] ,  (26) 

which  resulted  in  Eq.  (17)  of  this  work.  Similarly,  when 
there  is  a  driving  electrode  and  the  composite  plate  is  driven 
by  a  voltage  Ve" ",  from  Eqs.  (3.26)  of  Ref.  7  the  solution 
takes  the  form 

u[  =  u(  -  (e{,  Fx,/c,',/t  ')<>'"', 

if  1  =  ieU/cU  >«,'  +  ICx,  +  K  +  (  Vxx/h  ')  ] e  ",  (27) 

where  u{  is  given  in  (24)  and  u\  is  still  given  by  (25).  This 
solution  resulted  in  Eq.  (21 )  of  this  work. 


IV.  STRAIGHT-CRESTED  DISPERSION  RELATIONS 

In  this  section  we  obtain  the  straight-crested  dispersion 
relations  for  the  composite  plate,  the  thin  diaphragm  with¬ 
out  the  film,  and  the  bulk  semiconductor  because  these  de¬ 
termine  the  solution  functions  that  are  used  in  the  variation¬ 
al  condition  (16).  Since  in  this  work  we  are  considering  strip 
electrodes  as  shown  in  Figs.  1  and  2,  we  are  interested  in  the 
straight-crested  eigensolutions  varying  with  x,  in  each  re¬ 
gion.  In  the  composite  region  we  specifically  treat  propagat¬ 
ing  waves  when  the  ground  electrode  is  and  the  driving  elec¬ 
trode  is  not  present  because  these  are  somewhat  simpler  ( less 
cumbersome)  to  treat  than  propagating  waves  between 
shorted  electrodes,  and  on  account  of  the  small  piezoelectric 
coupling  and  the  small  wasenumbers  of  interest  in  the  com¬ 
posite  region  of  the  resonator,  the  dispersion  relations  and 
the  associated  solution  for  the  case  of  shorted  electrodes  can 
readily  be  obtained  from  those  for  the  case  considered  here 
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simply  by  employing  the  appropriate  relations  in  Sec.  Ill  for 
the  pure  thickness  solution  in  both  cases  exactly  as  in  Ref.  7. 

From  the  differential  equations  and  boundary  condi¬ 
tions  presented  in  Sec.  II  we  see  that  for  the  straight-crested 
waves  varying  along  x , ,  which  are  of  interest  in  this  work,  we 
maytakeu,  =  Oandfrom  Eqs.  (1),  (7), and  (5)  wehavethe 
nontrivial  differential  equations 


c(\  u\.u 

+  (fj\  4  C474  )«{,., 

—  PJ“{' 

(28) 

C4XII 

+  (fp  4-  c4'4  )«/,, 

+  *3.1  M/.33 

+  eL<f 

.11  =/ 

0  'u 

(29) 

-  /, ,  =  0, 

(30) 

fii«l.ll 

+  ( Cl ,  +  C44  )u\ ,  , 

C44  w  1 33 

=  P'u\, 

(31) 

^44  ^  3.M 

+  (<*',«  + 

+  <1.1 

(32) 

and  from  Eqs.  (9)-(  11),  ( 

3 ) ,  ( 6 ) ,  and  ( 8 ) ,  we  see 

that  the 

boundary  conditions  for 

the  unelectroded 

plate 

with  a 

ground 

electrode  take  the 

form 

+  c{,u( ,  4-  e(x(f 

=  0  at  .x ,  = 

-  h  ', 

(33) 

w/..i  +  i 

W x . j  —0,  c[xu{x  - 

-  f  - 

0  at  Jr , 

-  h  '. 

(34) 

fu«l,l 

+  c(,m{,  + 

-  r'1(tv', 

,  = 

"p"u\. 

IV  {  =  1/ 

\  at  x ,  =  0, 

(35) 

<44  («i' 

,  + «(,)  -ci,(u\ 

1  +  u'vi  1  = 

h  "p"u 

1  • 

u(  —  u 

1, 

'  =  0  at  .x , 

-  0. 

(36) 

4-  =  0. 

M|,i  + 

tv', ,  =  0  at  .x ,  =  - 

h\ 

(37) 

As  in  the  case  of  pure  thickness  vibrations,  in  order  to 
satisfy  Eqs.  (30).  (34),.  and  (36)„  we  first  take  if  'in  the 
form  given  in  Eq.  (26),  which  we  note  also  satisfies  the  con¬ 
dition  if  '  =  0  at  x,  =  0,  and  substitute  into  Eqs.  (29),  ( 33 ), 
and  ( 35) ,  to  obtain 


Cw  «  ui 

4  (£•,', 

+  ^44 

+• 

"  P  'u 

(38) 

+  c(xu( 

,  •—  0  at  ar , 

=  h 

(39) 

‘’i'.w/i 

-  f’u«Vl  = 

/t  >'m\ 

at  ,x ,  =  0. 

(40) 

As  a  solution  of  Eqs.  (28),  ( 38 ),  ( 31  ),  and  ( 32 ),  we  take 


-  (Ay- 

'  +  *>  )e'"r'"'. 

(41 ) 

»  )(FV', 

(42) 

in  which  the  subscripts  a. h  take  the  values  1  and  3  but  skip  2. 
The  solution  functions  in  (41)  satisfy  Eqs.  (28)  and  (38) 
provided 

ir(tA  +  f A  l  0.  o A  (  f  a[,  A  {  ■-  0,  (43) 
t? !,  B  (  <r{,B  l  0.  rr',  B  I  f  < j'„B'x  0.  (44) 


where 

a'\\  =  “  Cm  5  2  -  <■«  V]  +  P  V, 

=  -  (fl)  +  CL  )§T)/’ 

a{>  =  -  C44S  vj  +pJ<o2.  (45) 

and  the  solution  functions  in  (42)  satisfy  Eqs.  (31 )  and  (32) 
provided 


(7  i  |  A  , 

+  <7 

1  <  A  , 

=  0, 

0  ,  A  \  4-  <7  V. 

A  \ 

=  0, 

(46) 

- 

—  <7'l 

iB\  = 

=  0, 

—  a  ,  B  |  +  £7 

,B\ 

=  0, 

(47) 

where 

^  n 

=  ■ 

-  <7'i  1  i 

2  ~  cl,  T);  +  p'o)2. 

<7',, 

=  - 

-  (<7'n 

+  C14 

)i7,. 

<7,', 

=  - 

-  C44S  : 

'  -  c\ 

1  v;  +  p'w2. 

(48) 

Each  of  these  four  systems,  Eqs.  (43).  (44),  (46). and  (47). 
of  two  linear,  homogeneous  algebraic  equations  in  two  am¬ 
plitudes  yields  nontrivial  solutions  when  the  determinant  of 
the  coefficients  of  the  amplitudes  vanishes.  Both  determi¬ 
nants  for  the  film  are  identical,  as  are  both  determinants  for 
the  layer.  Each  of  the  two  independent  determinants  is  qua¬ 
dratic  in  i  :.  if,  and  m:.  Hence,  for  a  given  £  and  to,  each 
determinant  yields  two  independent  r/(  t;J 1  \i)'f  and  1/' 1 
»/::i)  and  each  t/‘ ' ’  yields  independent  amplitude  ratios 
from  either  of  the  two  equations  leading  to  each  of  the  four 
determinants.  Let  us  denote  the  eight  sets  of  amplitude  ratios 
by 

A  !'  =  P  "A  B  C  =  x"B  A  \'  =  p‘A  v, 

B  ]'  =  v  B  V,  /  =  1,2.  (49) 

As  a  solution  of  the  boundary  conditions  (  34),,  (39), 
(40),  (35),,  (36),  ,,  and  (37),  we  take 

u{  =  £  (A  +  B  [me  . , 


u{  =  £  (//  ’’"A  4-  \  "nB  {me 

•n  I 

u\  —  £  (.4  +  B\me  ">■■■■' 


u\  =  £  (p'"'A  're'"  -'  +  v-B\"'c  (50) 

•>,  1 

Substituting  from  Eqs.  (50)  into  the  boundary  conditions 
(34),.  (39),  (40).  (35),.  (36),.,,  and  (37)  and  employing 
Eqs.  (49).  for  h  "  =  0  we  obtain 


2242 


J  Appl  Phys  .  Vol  60.  No  7,  1  October  1986 


Shick,  Stevens,  and  Tiersten 


2242 


I  u 

n  l 

1  +  c 

:<  )e"‘ 

,.,h  ’ 

-t- 

)e 

"\<nB  ' 

]=0, 

i  u 

n  I 

4 

)e”'"  h  '  w 

B 

;■( 

~  W*  +  $)e  '""h 

]  = 

0, 

i  \a 

n  i 

fn ,  „  '  »-  In 

.  <P|.s./4 

4  c 

'■  >  + 

B 

"(c 

'f'T 

1 

1 

”'(c. 

1  > 5/' +  c\ ,  if  )  -  B'"( c; ,  o-'"  -  c\ ,  j;,,,  )  ] 

v  (a  ; 

■'  4  b  r  - 

a  \- 

-  b  r ) 

=-  0 

II 

V  [f 

1 

"<\a  1  V-../4 

4  >  *  b  : 

)( 

-  4  i  )  A  \V 

44  <  '/ 

P  ~  B  \"c'M  (  -  )/  ,,  -4  c  )  ]  --  0. 

V  <//  ' 

a  r  -  >  "■ 

B 

B  "' A  I" 

r" 

b  ■;•)  =0, 

i  !  < . 

»  c,. 

)c 

V  .  *• 

4- 

B  \"(c -  <•;,  7/,„ 

)  e  ''  ’ 

-  0. 

I  [ » 

'  1  >1  .ft  '  ^ 

2')c 

•  v  .  •> 

B  . 

'( 

4  C)e"  -h  ] 

0. 

Equations  (51)  const i lull.'  a  system  of  eight  linear  homogen¬ 
eous  algebraic  equations  in  A  B  ",  A  and  B  ,,  which 
>  leld  nontrivial  solutions  when  the  determinant  of  the  coeffi¬ 
cients  vanishes  For  a  given  geometry  and  <■>  the  resulting 
equation  yields  an  infinite  number  of  roots  i\  each  of  which 
determines  a  point  on  the  dispersion  spectrum  and  yields 
amplitude  ratios  front  any  seven  of  the  eight  equations  in 
Fqs.  (51  )  (  alculations  are  performed  by  first  selecting  val¬ 
ues  of  o  and  c.  which  enable  the  determination  of  the  t/,  and 
V  i  from  the  two  independent  2  •  2  determinants  obtained 
front  F-.qs  ( 42  )  or  ( 44 )  and  Fqs.  (46)  or  (47 ),  respectively. 
Iltese  values  of  >/,  and  i/„  enable  the  determination  of  the 
amplitude  ratios /j  \  v  and  v"  from  either  of  the  equa¬ 
tions  in  Fqs.  (42),  (44),  (46),  and  (47),  respectively  Then 
all  quantities  in  the  determinant  obtained  from  Eqs.  (  51 )  are 
known  and  the  resulting  dcterminantal  equation  either  is  or 
is  not  satisfied.  If  it  is  satisfied  the  v  alues  ofru  and  c  selected 
are  correct  and  constitute  a  point  on  the  dispersion  curves.  If 
not,  change  either  «>  or*  and  repeat  the  calculation  until  the 
boundary  condition  determinantal  equation  is  satisfied. 

I  he  pertinent  dispersion  curves  for  the  aluminum-ni¬ 
tride  film  on  the  gallium-arsenide  layer  composite  plate 
have  been  calculated  and  are  shown  in  Fig.  4.  We  do  not 
distinguish  between  the  electroded  film  with  shorted  elec¬ 
trodes  and  the  uneleetroded  film  in  the  figure  This  is  done  in 
great  detail  in  Ref.  7  Furthermore,  when  required,  the  dif¬ 
ference  between  the  two  can  be  calculated  from  («>'  -  <7/  ) 
given  in  Sec.  Ill  The  dispersion  curves  are  for  a  film  thick¬ 
ness  of  7  //m  and  a  diaphragm  thickness  of  14 /tm 

In  the  region  of  the  bulk  semiconductor  the  solution  is 
considerably  simpler  than  in  the  composite  region  of  the  res¬ 
onator  because  there  is  only  one  section  with  traction-free 
upper  and  lower  surfaces.  However,  although  this  problem 
can  be  simplified  even  further  by  placing  the  coordinate  sys¬ 
tem  in  the  center  of  the  plate,  it  is  not  convenient  for  us  to  do 
tins  in  this  work  because  the  solutions  in  each  region  are  put 
together  in  the  variational  equation  (  16)  when  the  solution 
to  the  forced  vibration  problem  is  obtained  in  the  next  sec¬ 
tion  In  the  region  of  the  bulk  semiconductor  the  differential 


I - 

equations  are  identical  with  those  given  in  (21 )  and  (22), 
but  with  carets  over  the  variables,  and  which  for  clarity  we 
write  here  in  the  form 

('ii"i.n  4  (fjv  +  c44  )«;.,,  .  c\4u\  ,,  —f)'u\, 

4'  f^ii  4  cii  Iw’i.i <  •+•  c\> u',  ,,  =  p'u\.  (52) 

From  (6)  and  (15)  with  u,  —  0  the  boundary  conditions 
take  the  form 

=°  at  ar,  -  0  and  .tr,  =  -  it\  (52) 


I  l<>  4  Dispersion  curves  for  ihe  composite  region  consisting  of  an  alumi¬ 
num  nit  rule  lilm  on  a  gallium  arsenide  layer. 
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L"vv^ir«jr"ir»ywjnui'j(  'jrjrw/v  vw’J  ’.■v’ 


^44  («',,+  u',1  )=  0  at  x.,  =  0  and  x,  =  —  A'.  (54) 

As  a  solution  of  (52)  consider 

u,',  =  (A  :,cn'"''  +  B'ae  (55) 

which  satisfies  (52)  provided 

cr , , ^4  =0,  cr  ',,.4  ,  +  a  \  \A  \  =  0,  (56) 

—  (x\yB\  =  0,  —  er  -I-  a  \yB =  0,  (57) 

where  a  J,  ,  tr  and  a  are  given  in  (48).  Each  of  the  two 
systems  of  linear  homogeneous  algebraic  equations  in  two 
amplitudes  yields  nontrivial  solutions  when  the  determinant 
of  the  coefficients  of  the  amplitudes  vanishes.  Both  determi¬ 
nants  are  identical  and,  in  fact,  are  the  same  as  the  earlier 
determinants  for  the  layer.  Thus,  for  a  given  $  and  to  each 
determinant  yields  two  independent  r/,  (.!]{' )  and  each 
))["  yields  amplitude  ratios  from  each  system  of  linear  alge¬ 
braic  equations,  i.e.,  (56)  and  (57).  Hence,  there  are  four 
amplitude  ratios,  which  we  denote  by 

A  ,  r-fl  ‘A  V,  ft =--  v‘B  v.  (58) 

in  which  the//1'  and  r'‘  are  the  same  as  in  (49). 

Asa  solution  of  the  boundary  conditions  ( 53)  and  (54), 
we  take 

u\  t  ( A  VV”  a-  BT'v  )f'~ '  e""‘. 


u7/2ri(MHz) 


100  50 


(10s  m’1) 


100  150 


F  IG.  5.  Dispersion  curses  for  the  gallium-arsenide  laver  without  the  thin 
film 


«;  V  i(t'mA  re"-'-  4  v'mBTe  (59) 

I 

Substituting  from  (59)  into  the  boundary  conditions  (53) 
and  (  54  ).  we  obtain 


v  f  A  ) 

ft  \ 

•  « ',"(£•;. ftv"  <•'„>/,,,)  I  o. 


tions  are  performed  as  in  the  case  of  the  composite  region 
with  the  given  attendant  reduction  resulting  from  the  elimi¬ 
nation  of  the  film. 

The  pertinent  dispersion  curves  for  a  I4-//m-thick  gal¬ 
lium-arsenide  diaphragm  have  been  calculated  and  are 
shown  in  Fig.  5.  The  lowest  17  real  dispersion  curves  for  the 
bulk  gallium-arsenide  plate  have  been  calculated  and  are 
shown  in  dimensionless  form  in  Fig.  6.  For  the  fundamental 


X  M  -t  P  +•  ft  Y’(  til]  -  0, 


X  [  A  ( c, ,  t// t  f'u  7,4  )c  '' 

f l  I 

t  cr  "  -  Cm  7„,  "  ]  •-  0, 

X  [  .4  V(  7 ‘  c  )<’ 

i  #'/'<  7,,*'”  *  i' 1  0.  (60) 

Equations  ( 60)  constitute  a  system  of  four  linear  homogen¬ 
eous  algebraic  equations  in  I  .  and  H  .which  yields  nontri¬ 
vial  solutions  when  the  detei mutant  of  the  coefficients  van¬ 
ishes  At  this  point  it  should  be  noted  that  if  4  7  and  It  "  a  re¬ 
written  as  complex  conjugates,  the  sy  stem  can  be  simplifier! 
into  solutions  respectively  sy m mot r i  and  antisymmetric 
about  the  centerline  ol  tl'.,  plate,  which  may  be  treated  sepa- 
i  it  el  y  However,  tins  is  riot  pa.  ocularly  convenient  tor  us 
because  we  need  both  the  sv mmetric  and  antisymmetric  so¬ 
lutions  and  vve  have  already  treated  the  composite  plate, 
which  programs  are  i  caddy  modified  for  this  case  Clearly  , 
this  solution  holds  for  the  region  of  the  thin  diaphragm  with¬ 
out  the  film  provided  only  that  h  is  replaced  by  h  .  Calcula- 


0  ?0  4C  60  80 


(10*  nC'l 

I  -  It  >  h  Dispersion  curses  I'm  the  ^.il  hum  arsenide  wafer  with  It  Ihcduncn 
sinnlcss  (requeues  noiinuh/cU  ssith  res|s.-et  lo  the  tirsi  thickness-shear  fre- 
queues 
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essentially  thickness-extensional  resonance  for  the  compos¬ 
ite  dimensions  mentioned  (around  132  MHz),  this  number 
of  dispersion  curves  is  for  a  gallium-arsenide  wafer  thick¬ 
ness  of  about  5  mils.  In  this  work  we  perform  calculations  for 


wafer  thicknesses  up  to  8  mils  for  which  there  are  30  real 
dispersion  curves  for  a  frequency  of  132  MHz.  However,  we 
do  not  bother  to  show  the  figure  for  more  than  17  in  this 
work. 


V.  FORCED  VIBRATIONS  OF  COMPOSITE  RESONATOR 

In  this  section  we  determine  the  essentially  thickness-extensional  vibrations  driven  by  the  application  of  a  steady-state 
driving  voltage  to  the  strip  electrodes  on  the  surfaces  of  the  piezoelectric  film  of  the  composite  resonator  shown  in  Figs.  1  and 
2.  Since  we  include  radiation  into  the  bulk  semiconductor,  we  can  use  the  solution  to  calculate  the  Q  at  resonance  of  the  mode 
resulting  from  radiation  into  the  bulk  semiconductor. 

In  accordance  with  the  earlier  discussion  we  take  the  approximate  solution  in  the  form 

u'=  x  ft'*'*™'-  ^~i  sai,  a:  =  x  a" 


fi  _  i 


H~\  C{yh 

■f-  X  u\  =  X  K,n,u'J"\  u'a= 


Uf,  =  2,  U'a  =  2-  *  K=  2.  (61) 

av  I  „  I  )-  =  1 

where  h,  n,  and  h  denote  the  number  of  branches  of  the  dispersion  curves  included  in  the  electroded  composite,  unelectroded 
composite  or  diaphragm  and  bulk  semiconductor,  respectively,  which  are  given  in  Figs.  4—6.  The  eigensolution  functions 
u'aP\  u'J"',  u“r>  denote  the  solution  functions  that  satisfy  the  homogeneous  differential  equations  and  the 

boundary  conditions  on  the  major  surfaces  for  each  of  the  respective  regions  and  which  were  presented  in  Sec.  IV  and  yield  the 
dispersion  curves  shown  in  Figs.  4-6.  At  this  point  we  note  that  since  the  driving  voltage  V  is  applied  over  the  region 

—  /  <  jc ,  <  /,  the  eigensolution  functions  in  the  electroded  region  of  the  composite  resonator  will  be  symmetric  about  .v,  =  0 
and,  consequently,  the  factor  in  (50)  will  be  replaced  by  (<?'*''  +  e  ‘-'])/2  in  the  expressions  for  u{  and  u\  and  (e'i<( 

—  e  '/2i  in  the  expressions  for  u(  and  u) ,  respectively,  in  (61 ), , .  Since  the  solution  functions  u'J'\  u u'J"\  and 
u'J r  1  satisfy  the  homogeneous  differential  equations  and  boundary  conditions  on  the  major  surfaces  for  each  of  the  respective 
regions  and  the  solution  is  symmetric  about  xt  =  0,  when  (61 )  will  be  substituted  in  ( 16)  all  that  will  remain  may  be  written  in 
the  appropriate  form 


') 


dx , 


nh  ’  cj/  Vx  rh>  -  — 

pW-~4e^6u{dx,dxl-  ( T {,  8u{  +  T {\  8u{ 

'  c{xh  *  Jo 

dx,+  — j  [ (r;,  -  r;,  )(Su\  +  Su\ 

>,  ,/  2  j  h  • 


+ 


J 


(T Su)  +  T) j,  Su\ ) 


f  (  T  | ,  —  /  ■, ; )  (Su\  4-  Su ,)  +  («!  —  u  (  )  (ST'n  +  <57^,  )  +  (u\  —  u, )  (6T\  ,  +  ST  u)],  /  dx , 

1 


4 


j  [  ( T  l\  -  T(\  )(8u{  +  8u{ )  4  (T/,  -  T{,)(Su{  +  Suf) 

4  (u{  -  u{)(8T{\  +  ST  (\ )  +  («{-  u{)(8T> ,  +  8T{, )  ] ,  ,  dx, 

4  ■— ■ ■  f*  [  ( 7  ']  j  —  T 1 1  )  {Su)  -f  <5i?|  )  4-  ( J  )  (  —  T\ ,)  {Su j  4  Su', ) 

2  Jh'  h' 

4  ( u  j  —  fi  j  )  ( <5  T  j  j  -f-  S  T 1 1 )  4  ( u ,  —  u,)  (ST) | !  +  <5F|,)]l  rf  dx,  =  0, 


(62) 


for  the  case  when  trapping  is  not  present  and  the  film  extends  over  the  entire  diaphragm,  as  shown  in  Fig.  1.  When  the  film 
ends  at  the  edges  of  the  electrodes,  as  shown  in  Fig.  2,  the  fifth  integral  in  (62)  is  not  present  and  the  second  integral  is  replaced 
by 


f 


(T{,6u{  +  T{]6u{)\Xi  j  dx 


(63) 


When  trapping  is  present  with  a  notch  under  the  electroded  region,  as  shown  in  Fig.  7,  the  fourth  integral  is  taken  from  —  h ! 
to  0  ( instead  of  -  !i  to  0)  and  there  is  an  additional  term. 


/; 


(  /  ;,  Su)  +  Su\ ):  ,  dx„ 


(64) 


in  Eq.  (62 ). 

Since  the  solution  functions  in  (6) )  are  fixed,  only  the  A' l/l',  K  and  K  ' y>  are  varied  when  (61 )  is  substituted  into  (62). 
Accordingly,  substituting  from  (61 )  into  (62),  employing  (3),  (6),  and  (50)  with  the  appropriate  aforementioned  replace¬ 
ment  in  (61 ),  ,  and  (56),  and  performing  the  integrations,  we  obtain  an  equation  of  the  form 


2245 


J  Appl  Phys  .  Vol  60,  No  7.  1  October  1986 


Shick,  Stevens,  and  Tiersten 


2245 


y<5*‘*>^>  +  y  y  6K WK "%<■>  +  X  X  <sa:("!a:uv  +  y  y 

0  —  1  0=15-1  (1  -  I  <  -  I  y  ^  1  £  1 


+  X  X  iSK^K^blla  +5K^K^b0„)+  y  y  (6K'a'k,r>bur  +  8Ku,Kla,b„r)  =0 

0=  I  o  =  1  a  -  1  y  =  I 

_  A 

where  the  expressions  fora/M,  a„f,aH,  bPn,bliy,  and  ciP)  are  given  by 

=  {  I  I  f  VfDo'hfLj  +g(v;,,-:nPrX’h%n.i  +  g(tf,  ~  V'j'X’hfX 

+  *{  ~  V“m  -  V%  )n  'h  ",La  +  g(  v':„,  +  rfl  )"  b  fn +  glr£m  -  rf.„ )°  hM  f  , 

+  *<%  -  v“m  )".  /, .A f +  g(  -  jj"t  -  )°  *.A £,4  ] eH“ ' 

v.  =  ~  X  X  Ux;,,,  +  < ),* 'A X ,  +«(v;„  -  < >5  'A +  g(<  -  >7;;, >15  'h%n , 

/»?  1  /r  I 

+  g(  ■-  y:„  -  < ).'!  'a ;1.4  +  g(f/" ,  +  »/('„ )"  „ .a  ,  +  g(F/7m  -  1/7  )n  ,/*  :^„2 

+  *«  -  v'/.„ )"  ,  A +  g(  -  f/:„,  -  >/' )"  „  a ;;C„ 4  ] e',L  ‘ 

1  X  X  {U(  v;.„  +  y',1, ) 'a  ,  +  g ( v;:„  -  < ) 'a ;i ,  +  g i <  -  r/;„ ,  j *  h  %„ 

>t  i  «  i 


0/<„  = 


2  rn  1  „  I 

+  g(  -  Vj„ )!',  'A  ;;;:„4  ]e";>  •  ■- 

+  [s(v'7,„  +  y', )"  „  a;'1,  +g(r,''„  ->/('„)"  ,  /) +g(r/,„  -  )°  irh:X 

+  g(  -  y',,  -  »/'„ )"  *  a  i‘,,4  ]  * 4,,,/  +  «"*■■  ■ 

-  2  [g(  »/>«  +  v'l '( <;a  ft  +  ft )  4-  gt  7;;,:,,  -  ///,  >* '( <'0 ft  +  r (:;'b  ft ) 
+  g(<  -  >/7„  )o  Vft'a  ft  +  <r{”A  ft )+g(  -  v;„,  -  »/;, ),';  '«•/? -  a/;;^ /;; 


V-  =  —  X  X  [s(v'L  + *.A;t,  +g('/;:„  -  >/'„)" ,,  * 

^  m  I  n  I 


^C»  /,  a// 


+  g(v'L  -  v'D 


,0 


h« -  •=  y  X  X  {[*(»/"-  +  #»)"  /,  +^(v;:„  -  v",, )"  „,Aii2  +  gdfL  -  vzx  „m±, 

£  m  I  I  ' 


+  *<  -  Vl  vC, )"  „  a;14]  +2[g(7/;;„  <  ?y;,) 

+  g ( v::„  -  »y: )  ;;  ( n: r  ■;<  +  f*,?; ^  7;,  )+g(»y:  -  v;:„ )  •;'(</  v;; .4  +  *  7:; .4 ) 
+  -  v"«  -  >/!, )  ;;  <a-  7:;«  7/,  +  a  zb  z )  ] }, 

-  T^V’-777  7T-  X  -v£>] 


1 

,-tJ 

h,  . _ I 

r  ■  t  1  -f 

L  /  -4 

[. - 73 _ _ 

_ i 

V 

t _ 

In  (66)  and  (67)  we  have  employed  the  definitions 
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the  elect  r.  Mies  to  cause  trapping  m  the  fundimcntal  niotle  when  it  does  not 
trap  in  the  flat-plate  configuration 
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and  ?V",  r”".  <7  V" ,  <7  are  obtained  by  replacing  /  with  v  in 
( 70 ) .  The  expressions  for  the  unelect  roded  composite  or  dia¬ 
phragm  and  bulk  semiconductor  regions,  respectively,  i.e., 

t{"\  t{”\  cr£.  a(:\  r'Z .  r'Z .  a  Z  -  °  and  r'Z, ,  r\”l  ,  <7  Z,  ■ 

a  |”'  are  obtained  by  replacing  the  quantities  with  the  bars  in 
( 70)  by  the  associated  quantities  without  bars  for  the  unelec  - 
troded  composite  or  diaphragm  or  with  carets  for  the  bulk 
semiconductor.  Furthermore,  in  (66)  we  have  employed  the 
additional  definitions 
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(71) 

and  the  h ’s  for  the  other  regions  are  obtained  by  replacing  the 
quantities  in  (71 )  by  the  appropriate  quantities  associated 
with  each  region.  Since  the  variations  in  (65)  are  arbitrary, 
we  obtain  the  inhomogeneous  linear  algebraic  equations  for 
the  A' 


1  ,  A ,  and  K' 

1  in  the  form 

'  a. .  v  V  A 

'  =  -  I  F1' 

”  1 

a 

‘  •  V  A  ' 

a  1 . n. 


V  A  -  a  .  •  V  A  "h„,  —  0.  y  =--  1 . n.  (72) 

^  \  ■<  \ 

which  constitute  n  •  n  •  h  inhomogeneous  linear  algebraic 
equations  for  Ft  A  ,  hA"  .  andi/A'"1. 


VI.  QUALITY  FACTOR  RESULTING  FROM  RADIATION 

In  this  section  we  calculate  the  quality  (actor  due  to 
radiation  into  the  bulk  semiconductor  for  each  of  the  three 
eases  from  the  analysis  presented  in  Sec.  V.  Although  the 
solution  may  be  obtained  from  the  analysis  in  Sec.  V  at  any 
driving  frequency  for  which  all  the  pertinent  dispersion 
curves  are  available,  in  this  work  we  are  interested  in  the 
solution  only  at  t lie  fundamental  essentially  thickness-exten- 
sional  resonance  for  each  ol  the  three  configurations.  It  is 
clear  from  experience 1  1  that  under  these  circumstances  we 
need  consider  only  the  essentially  thickness-extensional 
branch  in  the  composite  region  ( either  electroded  or  not )  of 
the  resonator,  which  is  the  curve  labeled  1  in  Fig.  4.  How'- 
ever.  in  the  thin  region  of  the  semiconductor  without  the  film 
and  the  bulk  semiconductor  all  pertinent  dispersion  curves 
shown  in  the  respective  Figs.  5  and  6  must  be  included  to 
obtain  accuracy.  The  use  of  only  one  branch  for  the  compos¬ 
ite  region  means  that  we  always  have  n  -  1  and  n  —  1  for  the 
configuration  shown  in  F  igs.  1  and  7  but  not  for  the  configu¬ 
ration  shown  m  Fig.  2.  for  which  all  curves  shown  in  Fig.  5 
must  be  included. 

As  usual,  the  quality  factor  Q  is  defined  by 
Q  (A  *  U)/E"\  (73) 

where 
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1 1  k 


30CG* 
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MG.  X  Quality  factor  vs  wafer  thickness  when  trapping  is  not  present  lot 
the  composite  resonator  configuration  show  n  in  Fig  I 
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E"‘-  =  -2  frff-J-Re  f  (-r.Jj).,  „dx{.  (74) 

Jo  2  J  h 

in  which  T  is  the  period  of  the  v  ibration  and  we  hav  e  taken 
the  liberty  of  writing  the  integrals  in  (74) ,  *  over  discontin¬ 
uous  functions  to  achieve  brevity.  For  a  given  geometry  and 
mode,  resonance  is  determined  by  obtaining  Q  over  a  ( nar¬ 
row)  range  of  frequencies  and  finding  the  frequency  for 
which  Q  is  a  maximum. 

In  performing  the  calculations  we  have  found  that  it  is 
imperative  that  all  radiating  plate  waves  in  the  thick  region 
of  the  gallium-arsenide  be  included  in  order  to  achieve  accu¬ 
racy  Since  at  a  given  frequency  the  number  of  radiating 
waves  in  a  plate  goes  up  significantly  w  ith  thickness,  we  hav  e 
considered  gallium-arsenide  wafers  no  thicker  than  0  mils  at 
a  frequency  around  1 32  MHz.  for  which  there  are  30  radiat¬ 
ing  plate  waves.  Results  are  presented  for  wafer  thicknesses 
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FIG  1 1  Quality  factor  vs  wafer  thickness  when  trapping  is  not  present  for 
the  composite  resonator  configuration  shown  in  Fig  2 
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FK »  10.  Quality  factor  \s  distance*  from  edyc  of  electrode  to  edge  of  dia¬ 
phragm  when  trapping  is  present  for  the  composite  resonator  configuration 
shown  in  1  ig  T  The  width  of  the  electrodes  2/  varies  from  5X0  to  420 //m. 

ranging  from  4  to  8  mils  because  tins  is  considered  to  be 
within  the  practical  range.  All  the  results  presented  are  for  a 
film  thickness  of  7  // m  and  a  diaphragm  thickness  of  14/tnt 
and  the  lateral  dimensions  of  each  configuration  were  ad¬ 
justed  slightly  to  maintain  the  same  resonant  frequency  for 
computational  convenience.  Since  the  calculated  Q  is  a  very 
rapidly  varying  function  of  the  wafer  thickness,  calculations 
had  to  he  performed  for  very  small  increments  in  thickness, 
i.e..  1  ji m,  in  order  to  get  all  the  peaks  and  valleys  in  the 
interval. 

In  the  absence  of  trapping,  for  the  ease  shown  in  Fig  1 
the  results  are  plotted  in  Fig.  8,  which  shows  the  aforemen¬ 
tioned  sharp  sat  iation  in  Q  with  wafer  thickness.  It  can  be 
seen  from  the  figure  that  the  highest  Q  obtained  is  about 
4750  and  the  lowest  is  about  10.  and  there  are  about  10  peaks 
for  thicknesses  between  4  and  8  mils,  file  highest  valleys 
have  Q's  of  about  700  and  2(HX).  respectively.  For  the  ease 
show  n  in  Fig.  2.  the  results  are  plotted  in  Fig.  9.  which  shows 
variations  in  Q  with  wafer  thickness  similar  to  Fig.  8.  but  in 
this  case  the  Q ’s  are  considerably  higher.  It  can  be  seen  from 
the  figure  that  the  highest  ^obtained  is  about  1  5  OOOand  the 
lowest  Q  is  about  200.  Calculations  were  performed  when 
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FIG.  12  Typical  thickness  displacement  along  the  surface  of  the  composite 
resonator  for  the  configuration  shown  in  Fig  2. 

trapping  was  induced  in  the  fundamental  mode  by  notching 
the  diaphragm  under  the  electrode  as  shown  in  Fig.  7.  As 
expected,  the  Q  due  to  radiation  can  be  made  as  high  as  we 
w  ish  simply  by  extending  the  lateral  dimensions  of  the  film 
and  diaphragm.  Figure  10  gives  Q  as  a  function  of  (d  —  /), 
i.e..  the  distance  from  the  edge  of  the  electrode  to  the  edge  of 
the  diaphragm.  It  can  be  seen  from  the  figure  that  Q  in¬ 
creases  very  rapidly  with  (d  —  /)  and  at  d  —  21,  Q  =  70  000 
for  this  geometry.  For  these  calculations  the  wafer  thickness 
was  6  mils.  Even  when  trapping  is  present  the  Q  is  a  varying 
function  of  w  afer  thickness,  but  the  range  of  the  variation  is 
less  than  1/ 10  of  that  w  hen  trapping  is  not  present,  w  hich  is 
not  of  interest  since  the  Q  with  trapping  is  so  much  higher 
than  the  Q  without  trapping.  In  interpreting  the  foregoing 
information  it  should  be  remembered  that  the  high  Q' s  cal¬ 
culated  should  be  higher  than  the  actual  Q's  because  the 
material  Q  and  the  Q  due  to  radiation  into  the  air  are  not 
included.  The  results  discussed  above  mean  that  in  order  to 
obtain  reasonably  high  Q  when  trapping  is  not  present,  for  a 
given  wafer  thickness  the  thicknesses  of  the  film  and  dia¬ 
phragm  must  be  very  precisely  selected. 

Since  all  conditions  are  natural  conditions  in  the  vari¬ 
ational  condition  in  (62).  the  approximate  solution  does  not 
expressly  match  the  mode  shape  at  the  interfaces  unless  it  is 
exact.  Consequently,  the  extent  to  which  the  calculated 
mode  shape  matches  at  the  interfaces  gives  an  indication  of 
the  accuracy  of  the  approximate  solution  obtained.  Typical 
plots  of  the  //.-displacement  field,  which  is  the  large  one  for 
the  essentially  thickness-extensional  modes  considered  here, 
for  the  configurations  shown  in  Figs.  1  and  2  are  shown  in 
Figs  1 1  and  12.  respectively.  It  can  be  seen  from  the  figures 
that  the  u  ,-displacement  field  matches  quite  well  at  the  inter¬ 
faces.  Consequently,  we  can  conclude  that  the  approximate 
solution  obtained  is  quite  accurate. 
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1.  Introduction 
I  2 

In  a  recenc  analysis  ’  of  the  piezoelectric 
thin  film  on  semiconductor  composite  resonacor 
vibrating  in  an  essentially  thickness-extensional 
mode  it  was  shown  that  the  quality  factor  (Q)  due 
to  radiation  into  the  bulk  semiconductor  is  a  very 
sharply  varying  function  of  the  ratio  of  the 
thickness  of  the  resonator  to  chat  of  the  film  if 
trapping  is  not  present,  but  not  if  trapping  is 
presenc.  The  treatment  employs  a  very  accurate 
but  extremely  cumbersome  variational  approximation 
technique  and  is  restricted  to  the  case  of  strip 
electrodes  and  diaphragms. 

In  this  work  a  perturbation  analysis  of  the  Q  due 
to  radiation  into  the  semiconductor  wafer  is  pre¬ 
sented,  which  is  considerably  less  cumbersome  to 
use  than  the  earlier  variational  treatment  and  is 
not  restricted  to  the  case  of  strip  electrodes  and 
diaphragms.  In  the  treatment  the  resonant  mode  of 
interest  is  determined  from  the  equation  for  trans¬ 
versely  varying  essentially  thickness-extensional 
inodes  in  composite  resonators^  and  simple  approxi¬ 
mate  but  very  accurate  conditions  at  the  edges  of 
the  diaphragm.  This  resonant  mode  is  then  used  to 
determine  the  near  field  radiation  into  the  semi¬ 
conductor  wafer  by  means  of  a  variational  approxi- 

it 

mation  procedure  .  Of  course,  as  in  the  earlier 
1  2 

work  ’  ,  all  radiating  plate  waves  in  the  wafer  are 
included  in  order  to  achieve  accuracy.  Finally, 
the  mode  in  the  composite  resonator  and  the  radia¬ 
tion  field  in  the  wafer  are  employed  in  a  perturba¬ 
tion  integral'*  to  calculate  the  Q.  Calculations 
are  performed  for  the  cases  of  rectangular  elec¬ 
trodes  and  diaphragms  and  strip  electrodes  and 
diaphragms  both  when  trapping  is  and  is  not 
present.  For  the  case  of  strip  electrodes  and 
diaphragns  the  calculated  results  are  shown  to  be 
in  good  agreement  with  the  earlier  more  cumbersome 
1  2 

direct  calculations  '  . 

2.  Perturbation  Procedure 

Since  the  coupling  is  small  in  the  piezoelectric 
thin  film,  we  need  consider  only  the  elastic  por¬ 
tion  of  the  equation  for  the  first  perturbation  of 
the  eigenfrequency,  which  may  be  written  in  the 
form"* 

V  *  VH .  ^  >  a) 


where  w  and  tu  are  the  unperturbed  and  perturbed 
eigenfrequencies,  respectively,  and 

«-•  Jni(Tij8ruj^)dS-  <2> 

where  denotes  the  outwardly  directed  unit  normal 
to  the  surface  S  of  the  resonator  and  the  normal' 
ized  displacement  field  is  defined  by 

gu-uu/N  ,  N2  -  f  pulTuudV,  (3) 

J  J  J  )  ]  ' 

V 

in  which  u1!  denotes  the  purely  real  mechanical 
displacement  field  associated  with  the  mode  of 

interest  in  the  composite  resonator  and  T“ .  denotes 

i  i 

the  associated  stress  tensor.  A  cross-section  of 
the  composite  resonator  attached  to  the  semicon¬ 
ductor  wafer  is  shown  in  Fig.l.  A  schematic  plan 
view  of  the  assumed  composite  resonator  is  shown 
in  Fig. 2  along  with  the  assumed  approximate  edge 
condition  (Uj*0).  Thus,  for  the  problem  at  hand 

the  surface  S  in  (2)  denotes  the  surface  along 
which  the  semiconductor  portion  of  the  composite 
resonator  abuts  the  wafer.  The  calculated  mode  in 
the  compos ite  resonator  shown  in  Fig. 2  results  in 
tractions  t^  and  displacements  u^  that  are  applied 

to  the  wafer  along  the  surface  at  which  it  abuts 
£he  resonator  as  shown  in  Fig. 3.  These  tractions 
t^  and  displacements  u^  then  cause  radiation  fields 

in  the  wafer,  which  are  calculated  by  means  of  a  ^ 
very  accurate  variational  approximation  procedure  . 
These  radiation  fields  in  the  wafer  produce  a 
reaction  back  on  the  composite  resonator  along  S. 

The  stresses  T.  .  and  displacements  u.  in  (2)  denote 
ij  1 

this  back  reaction  field. 

Since  the  radiating  fields  in  the  wafer  are  chosen 
to  satisfy  the  differential  equations  in  and 
boundary  conditions  on  the  major  surfaces  of  the 
wafer  exactly,  all  that  remains  of  the  appropriate 
variational  principle  in  which  all  conditions  are 

natural  conditions  is  given  by^1 

J  (Tz'nk\i)SSids  +  J  VVV^k/5*0-  (4) 
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fcv: 


In  (4)  c.  and  u,  are  known  from  the  resonant  eigen- 
l  L 

solution  T*t,  .  ut  and  u.  and  T.  .  are  found  from  the 
ik’  1  l  ki 


expressions  for  the  solution  field  radiating  into 
the  yafer.  At  this  point  we  note  that  the  reaction 


u  ,  T^£  from  the  radiating  field  is  complex,  the 


real  pare  of  which  yields  a  small  change  in  fre¬ 
quency  when  substituted  in  (2),  which  is  negligible 
and  not  of  interest  here,  and  the  imaginary  part  of 
which  yields  the  attenuation  due  to  radiation  into 
Che  wafer  from  (2)  which  is  sought  here.  Frcxn  the 
well-known  relation  ui  *  u)  -  ix  / 2Q  and  (1),  we  obtain 


1  -  ix  /  H 


(5) 


Transversely  Varying  Thickness-Excensional 
Modes 


As  noced  in  Sec. 2,  a  plan  view  for  the  determina¬ 
tion  of  the  mode  of  interesc  in  the  composite 
resonacor  is  shown  in  Fig. 2,  in  which  the  -  denotes 
the  electroded  region,  the  S,  the  side  region,  the 
T,  the  top  (and  botcom)  region  and  the  C,  the 
corner  regions.  This  notacion  is  essential  for  the 
treatment  of  Che  trapped  modes,  but  is  not  needed 
for  the  untrapped  modes  because  in  the  uncrapped 
case  the  edges  of  the  electrodes  are  relatively 
unimportant.  It  has  been  shown  in  Sec. 5  of  Ref. 3 
that  the  homogeneous  equacion  governing  the  essen¬ 
tially  th icknes s-extens iona 1  modes  may  be  written 
in  che  form 


"*  n  n  n 

Mn(tT  +  Tt) 

cx 


-f  ;2  ,n  tyn 
c  f  *  s  f  1 

33  fn 


(6) 


1 


where  n  denoces  the  order  of  the  pure  thickness 
mode  in  che  composite  resonacor.  both  the  super¬ 
script  and  subscript  f  denote  che  film  and 

in  an  electroded  regionwhile 


in 
in  an 


in  "  'fn  fn 

unelectroded  region.  Furthermore,  it  has  also  been 
shown  in  Sec.V  of  Ref.  3  chat  at  an  interface  between 
electroded  and  unelectroded  regions  we  have  che 


continuity  of  fn  and  dfn/dn,  where  d/dn  denoces  the 


normal  derivative.  Within  che  aoproxima Cion  madeJ 


6-8 


in  obtaining  Eq.(6)  it  may  readily  be  shown  that 
for  either  completely  free  or  completely  fixed 
conditions  (here  completely  fixed)  along  che  edges 
of  the  diaphra^n  shown  in  Fig.  2,  the  appropriate 

*n  -  "  The  expression  for  the  coeffi- 


condicion  is  f 


cient  M  ,  which  is  very  important  because  its  sign 


indicates  whecher  trapping  is  or  is  not  present  in 
the  flat  composite  plate  configuration,  is  given  in 
Eq.  (4.44)  along  with  (4.41),  (4.42),  (4.37),  (4.38) 
and  (4.33)  of  Ref.  3  and  clearly  is  much  too  cumber¬ 
some  to  present  here’. 


It  has  been  shown  in  Sec.V  of  Ref.  3  that  the  dom¬ 


inant  u^-displacemenc  field  accompanying  the  mode 


is  given  by 


fn  fn ,n  , 

J3  *oU3  f  (VX2' 


c), 


sn.n, 

ou3  f  (x1,x;,t), 


(7) 
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where  the  superscripts  f  and  s  denote  the  film 
and  semiconductor  diaphragm,  respectively,  and 


represents  the  thickness  dependence  only  and  is 
given  by 
fn 
oU3 

f 


o  3 


,  fin  o  ,  „fln  o 

'A,  cos  TL_x,  +B.j  sm  Tl£nX3 


'fn  3 


0  <  x^  <  h 


sn 

oU3 


, sin  o  _sln  o 

'A,  cos  T1  x,  +B,  sm  T  x,  , 
3  sn  3  3  sn  3 


0  >  x3  > • 


(8 


and  che  amplitudes  are  given  in  Eqs.(5.20)  of 
Ref. 3.  Since  only  the  semiconductor  diaphra^n 
abuts  the  wafer  as  shown  in  Fig.l  only  the  variable 
in  the  semiconductor  are  relevant  here.  However, 


in  order  to  calculate  che  tractions  t^  and  displace 


oents  u^  that  the  mode  in  the  composite  resonacor 


exerts  on  the  wafer  as  shown  in  Fig.  3,  we  need  the 


s  s 

associated  displacement  fields  u£  and  u^  in  che 
semiconductor  diaphragm  in  addition  to  u*.  These 


functions  differ  for  che  different  f  and  associ 


ated  u§  ,  which  occur  in  the  different  regions  for 
the  trapped  case  but  not  for  the  untrapped  case. 


Since  the  f  has  the  same  functional  form  over  the 
encire  composite  resonator  in  the  uncrapped  case  as 
it  has  over  the  electroded  region  in  the  trapped 

case,  we  write  u^  and  u^  for  this  functional  form 

n  3  1 

for  f  only  .  It  ha  s  been  shown  in  Sec.V  of  Ref. 3 
that  for 


f 


_  —  ixt 

cos  §x3  cos  vx,e 


(9) 


which  is  che  above  mentioned  fn 


(Xj)  sin  gx£  cos  vx,  e 


ixt 


VS.  —  -  ixt 

~  t  (x  )  cos  §x.  sm  vx e 

c 


(10) 


where 


s  .sl  „o  ,  so  ^„sl 

r  =A3  sin  T]sX3 +A£  sin  X  Tt_x,  +  B,  cos 


+  B£  sin  k 


Vs*  * 


s  3  1  s  "3 

—  I  —  "> 

r+v‘ , 


(in 


and  x  and  the  amplitudes  are  given  in  (4.30)  and 
(4.26),  respectively,  of  Ref. 3. 


As  already  noted  ,in  che  untrapped  case  for  che 
fundamental  mode  for  any  n,  fn  is  given  by 


—  _ 

f  “  cos  gx£  cos  vx,  e 


given 

ittt 


(12) 


where  from  the  edge  conditions  for  the  mode  of 
inceros  t 


P  -  n/2d 


n/  2w 


(13) 


with  which  the  unperturbed  resonant  frequency 'x 
8  “ 


may  be  obtained  from 


(14) 


(19) 


f  ->  £  -•>  —  •>  2 

«*  -J"C33r‘f  +  V>  +  * 

In  the  Crapped  case  Che  expressions  for  f11  in  each 
of  the  respeccive  regions  shown  in  Fig. 2  are  given 
8 

by  _  _  „  -5(x,-i) 

"f11  "  E  cos  ?x^  cos  vx,  ,  f  S  *  (E$e 

s  5(xr£) 

£+e  )  cos  vx, 

_  -v(x^-b)  v(x,-b) 

f  -  (E  e  ‘  +Ele  ‘  )  cos  §x.  , 


nC  C  *5<V1)-V(X 
f  -Eie 

c  ?<x1-*)-v(x,-b) 


,-b)  c  -^(Xj-X^+vfx^-b) 

+  z2e  ‘ 
r  £(x..-i)+v(x,-b) 

. 


»  0  at  *  0  and  -  h 

In  considering  waves  radiacing  in  Che  +x. -direction 

—  g  *  s  ^ 

we  first  noce  that  since  v  is  small  ,  T^,  is  an 
order  of  magnitude  smaller  chan  and  u* ,  and  u!j 

is  an  order  of  magnicude  smaller  chan  u^  and  u*. 

,  S  *  $ 

Hence,  u,,  Che  x^-di f ferencial  equacion  and  T  are 
negligible  for  radiation  in  the  x^-direccion.  Of 

course,  equivale’  r  statements  hold  for  radiaCion 
in  the  x, -direction  since  £  is  small.  Accordingly, 

for  waves  with  slowly  varying  crests  and  propagat¬ 
ing  in  Che  -me^-direcrion  in  Che  near  field,  we 
write 


where  £  and  v  are  determined  from  Che  lowest  roots 
of  Che  transcendental  equations 


C(m)(g+(m)e 


?  can  cX  »  § 


v  can  Vi  *  v 


(l4e-^(d-X)) 
,  *2£(d-£) 


n^e-2v(w-b)> 
,  -  2v(w-b ) 


with  Che  aid  o;  (6.7)  of  Ref.  3.  All  amplitude 
coefficients  in  (15)  are  known  in  terms  of  E  from 

the  relations'" 


E  cos  --i 
=  2-(d -*) 


E  cos  vb 
=2v(w-b ) 


C  _ E  cos  ij  cos  vb _ _ 

1  '  (1.e-2i(d-i))(1.e-2viw-o)) 


- 2v(w-b )  C  _C 

Er  E3 


e-2?(d-i)EC 


£C  .  e-2t(d-X)-2v(w-b)£C  _  (1, 

S  S  ”TJ 

As  already  noted,  u^  and  u^  are  known  from  f  in 

the  electroded  region  from  the  relations  (10)  and 
(11)  and  for  the  other  regions  equivalent  rela¬ 
tions  which  are  not  shown  here  are  employed  . 


„  ,  ,  -”v  'x,  -if (x, -d)  . 

.  t-  (m)  3,  s  1  lxt 

+  ca  e  >e  e  ,  a  -  1,3,  (20) 

for  either  the  untrapped  case  or  the  electroded 
region  for  the  trapped  case  and  we  do  not  bother 
to  write  the  solution  for  any  other  regions  for  the 


tapped  case.  The 


i(m)  g+tm)  e-(m) 


determined  by  satisfying  (18)  and  (19)  in  the 
8  2  3 

usual  way  ’  ’  .  From  this  solution  the  lowesc  17 
real  dispersion  curves  for  the  bulk  gallium- 
arsenide  plate  have  been  calculated  and  are  shown 
in  dimensionless  form  in  Fig. 4.  For  the  fundamental 
essentially  thickness-extensional  resonance  of  the 
composite  resonator  consisting  of  a  7  js  thick 
aluminum-nitride  film  on  a  14  yn  thick  gallium- 
arsenide  diaphragm,  which  is  around  132  MHr,  this 
number  of  dispersion  curves  is  for  a  gallium- 
arsenide  wafer  thickness  of  about  5  mils.  In  this 
work  we  perform  calculations  for  wafer  thicknesses 
up  to  8  mils  for  which  there  are  30  real  dispersion 
curves  for  a  frequency  of  132  MHt.  However,  we  do 
not  bother  to  show  the  figure  for  more  than  17. 

Similarly,  for  waves  with  slowly  varying  crests  in 
the  x^-direction  and  propagating  in  the+x^- 

direction  in  the  near  field,  we  have 


4.  Variable-Crested  Waves  in  Wafer 

In  this  section  we  obtain  the  solution  functions 
for  the  near  field  waves  with  slowly  varying  crests 
and  the  associated  dispersion  relations.  The 
stress  equations  of  motion  and  the  linear  elastic 
constitutive  equations  for  the  semiconductor  wafer 

.  .  *  8,10,11 

may  be  written  in  Che  censor  form 


is  s.s  cs  s  «s 

ij,i  *  uj  ’  ij  CijkfUk,i 


where  the  notation  is  conventional,  and  we  note 
that  the  equations  are  too  cianbersonte  to  write  out 
in  detail.  Prom  eicher  Figs.l  or  3  we  see  that 
the  boundary  conditions  on  the  major  surfaces  of 
the  semiconducting  plate  may  be  written  in  the 
form 


2  i-lnUx 

- s  7  r-  *<ra)  -+(o)  3 

Ub  *  C0S  5xl  i  C  (tb  e 

tn*l 

-  fml  ■iVB')X1  “is  (X,*W ) 

+  8^(  }e  3)e  ‘  ellut  ,  b-2,3,  (21) 

for  the  same  types  of  regions  as  in  the  previous 
case.  From  this  solution  we  obtain  the  propagating 
dispersion  curves  for  near  field  radiating  waves 
in  the  x, -direction.  The  dispersion  curves  are 

just  like  those  in  Fig. 4. 

5.  Variational  Analysis  of  Radiation  into  Wafer 
In  this  section  we  determine  the  waves  radiating 
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into  the  semiconductor  wafer  due  to  the  tractions 
and  displacements  resulting  from  the  mode  in  the 
composite  resonator  by  means  of  a  variational 
approximation  procedure,  as  indicated  in  Sec. 2. 
Since  Che  near  radiation  fields  emanate  from  the 
edges  of  the  diaphragm,  which  are  normal  to  x^  and 

x^,  respectively,  we  may  determine  the  near  fields 
radiating  in  the  x^-  and  x2'directions  separately. 

Accordingly,  for  radiation  in  the  x  -direction, 

8  1 

from  (4),  we  obcain 


J  -ll  <* 


S  .aS  ,  2LS 

USu1+T13 


-Gs 


dx3  ' 


In  a  similar  way  for  radiation  in  the  x,- direction, 
we  obtain 


N 

z 

Y-l 


.5. 


( 26 , 


Equations  (25)  and  (26)  each  constitute  N  inhomo¬ 
geneous  linear  algebraic  equations  for  the  N  ur.- 

t.  /  ^  m  *  (%  \ 

knowns  from  (25)  and  che  N  unknowns  K_ ' ' 

frca  (26) ,  respectively.  When  the  3rd  ir 

-  {y  )  s 

have  been  determined,  the  near  field  racia- 

g 

ting  solution  is  known  . 


l: 

i 


i 

« 

.Y 

.1 

.1 

j 


13 


13' 


,*S  1 1  , 

.  dx_ 
3  ’ ?  Xj,  *d  3 


-n 

o 


dxj. 


(22) 


_s  ,  s 
are  T,  ,  and  g 


where,  as  already  indicated,  the  forcing  terms 
resulting  from  the  mode  in  Che  composite  resonator 

ana  g^  .  For  radiation  in  the  x, -direction 

u«  simply  interchange  subscripts  1  and  2  and 
replace  w  by  d. 

For  radiation  in  the  x^-direccion  we  now  expand 

the  solution  in  the  wafer  as  a  sum  of  waves  with 
slowly  varying  crests  in  the  x, -direction  and 

propagating  in  che  Xj-direccion,  which  were  dis¬ 
cussed  in  Sec. 4  and  are  given  in  (20).  Thus 


-s  „  J  £ <»)-.<*> 


w  b 

o-l 


(22 ) 


where  each  of  the  uS^°^  are  of  the  form  given  in 

„  a 

(20),  N  denotes  Che  number  of  branches  of  the  dis- 

g 

persion  curves  required  for  the  N  propagating 


plate  waves. 
0  (cr) 


Since  Che  uS^CT^  are  fixed,  only  the 


K'“'  are  varied  when  (23)  is  substituted  into  (22). 
s 

Accordingly,  substituting  from  (23)  inco  (22), 
employing  (18).  and  performing  Che  integrations, 

8  ^ 

we  obcain  an  equation  of  the  form 


z  [i  «3  ^ 


r-o) . 


(24) 


3*1  a-*  l 


where  Che  expressions  for  A_^_  and  <2^  are  too 

lengthy  to  present  here.  Since  che  variations  in 
(24)  are  arbitrary,  we  obtain^  the  inhomogeneous 

*  (o) 

linear  algebraic  equations  for  the  K_  in  the  form 

ft 

*-Cg,  B  -  1 . :"l  .  (25) 

o-l 


Z 


6.  Quality  Factor  Resulting  from  Radiation 

In  this  section  we  calculate  che  Q  due  to  racia- 
tion  into  the  semiconductor  wafer  from  the 
analyses  presented  in  the  previous  sections. 
Accordingly ,  from  Secs. 3-5  we  see  chat  che  per- 

g 

curbation  integral  takes  the  form 
w 


x^»d 


(Vl3  * 


VMJ  [-Vl3+"ll^ 

-hs 

-Juh>x  .-d’dx3j  +  Jdxl[  J  :-=r23 

-hs 

+  ^22^2^x  ^ 3 fc 23  1  2^ ^ ^ 


Subscitucing  the  imaginary  parts  of  the  solutions 
determined  in  Sec. 5  into  (27),  we  obtain5 


\  I  .  ~  l  a 

-  [•m-N.  u_  )  i ,  , 

4_-  l  a  j  i.  j 

-hs  a-l 


H  *  2if  f  dx, 

U  L  J  2 
-w 

*  w  *  (or)-S  (cr)  ,_s  1  dx 

+  ^('NS  U1  >TlllxL-d  3 

*k !  x  [-s'',=r,>4 

-d  -hs  \ *  I 

.  w  r.OO-TCv )  s  1  .  1 

“  x.-w 

which  may  now  be  used  to  calculate  the  Q  from 

Of  course,  as  in  Che  earlier  work^’~  all  radia 
plate  waves  in  the  thick  region  of  the  gallium 
arsenide  are  included  to  achieve  accuracy.  Si 
at  a  given  frequency  the  number  of  radiating  w, 
in  a  plate  goes  up  significantly  with  chicknes 
when  trapping  is  not  present  we  have  considers 
gallium-arsenide  wafers  no  thicker  than  3  mils 
a  frequency  around  132  MHc,  for  which  there  ar 
30  radiating  place  waves.  Results  are  present 
for  wafer  thicknesses  ranging  from  4  mils  to 
3  mils  because  this  is  considered  co  be  within 
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practical  range.  For  the  untrapped  case  the 
results  presented  are  for  a  film  thickness  of 
7  microns  and  a  diaphragm  thickness  of  14  microns. 
Both  strip  diaphragms  600  microns  wide  and  square 
diaphragms  with  lateral  dimensions  of  600  microns 
x  600  microns  were  considered.  The  strips  are^ 
treated  for  comparison  with  the  earlier  work1’  . 

For  the  trapped  case  the  results  presented  are 
for  a  film  thickness  of  12  microns  and  a  diaphragm 
thickness  of  14  microns  because  the  trapping  is 
considerably  better  for  this  combination  of  thick- 

8 

nesses  than  the  other  at  the  second  thickness  mode  . 
For  chis  combination  of  thicknesses  the  second 
thickness  mode  is  around  250  MHz.  In  the  trapped 
case  both  strip  electrodes  500  microns  wide  and 
square  electrodes  with  lateral  dimensions  of 
500  microns  x  500  microns  were  considered.  The 
lateral  dimensions  of  both  Che  scrip  and  square 
diaphragms  was  varied  and  the  wafer  thickness  is 
6  mils . 

In  Che  absence  of  trapping  for  the  strip  case  Che 
values  of  Q  calculated  from  Eq.(5)  of  this  work 
are  plocted  as  the  solid  curve  in  Fig. 5,  in  which 
the  dotted  curve  from  Fig. 8  of  Ref. 2  is  also 
plocted  for  purposes  of  comparison.  It  can  be  seen 
from  the  figure  that  the  agreement  is  quice  good. 
However,  although  the  highest  Q's  calculated  in 
this  work  are  very  nearly  the  same  as  those 
obtained  in  the  earlier  more  cumbersome  direct 

calculation  ’  the  lowest  Q's  calculated  by 
means  of  the  perturbation  procedure  tend  to  be 
nearly  an  order  of  magnitude  higher  than  those 

calculated  by  Che  earlier  direct  procedure  We 

are  not  absolutely  sure  of  the  reason  for  this 
discrepancy,  but  there  are  two  possibilities.  The 
perturbation  procedure  might  be  tending  to  lose 
its  accuracy  for  low  Q  because  of  Che  increased 
radiation  or  the  resonant  frequency  might  not  have 
been  sufficiently  precisely  determined  by  means  of 

1  2 

the  earlier  direct  procedure  ’  for  the  accurate 
determination  of  che  lowest  Q  values.  However, 

Fig. 5  reveals  chat  the  highest  Q's  calculated  by 
the  perturbation  procedure  are  consistently 
slightly  higher  than  those  calculated  using  che 
earlier  direct  procedure,  which  tends  to  support 
the  second  possibility.  The  figure  also  shows 
that  the  location  of  the  peaks  and  valleys  of  Q 
with  wafer  thickness  determined  by  means  of  the 
perturbation  procedure  is  in  quite  good  agreement 
with  those  obtained  from  the  earlier  direct 
1  ? 

calculation  ’  .  Also  in  the  absence  of  trapping 
the  Q's  calculated  from  Eq.(5)  for  the  square 
diaphragm  are  plocted  in  Fig. 6.  It  can  be  seen 
from  the  figure  Chat  the  peaks  and  valleys  are  in 
essentially  the  same  positions  as  in  the  strip 
case,  but  that  the  Q's  are  considerably  lower, 
roughly  between  1/2  to  2/3  of  the  values  in  the 
strip  case.  This  is  as  expected  because  of  the 
radiation  in  two  orthogonal  directions  for  rect¬ 
angular  diaphragms.  When  trapping  is  present  the 
Q's  calculated  from  Eq.(5)  for  both  the  strip  and 
square  case  are  plotted  in  Fig. 7  as  a  function  of 
id  -  j£  > ,  i.e.,  the  distance  from  the  edge  of  the 
electrode  Co  the  edge  of  the  diaphragm.  It  can 
be  seen  from  che  figure  that,  as  expected,  Q 


increases  very  rapidly  with  (d  -  l)  and  for  the 
same  value  of  (d  -  i)  the  Q  is  about  twice  as  large 
in  che  strip  case  as  in  the  rectangular  case.  If 
a  film  thickness  of,  say,  8  microns  had  been 
employed,  the  required  (d  -  L)  for  good  Q  would  be 


much  larger 
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Figure  5.  Quality  Factor  Versus  Wafer  Tt.itrtr.tsi 
When  Trapping  Is  Not  Present  for  tne 
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Figure  6.  Quality  Factor  Versus  Wafer  Thickness 
When  Trapping  Is  Not  Present  for  the 
Rectangular  Composite  Resonator 
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Trapping  Is  Present 
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ABSTRACT 

In  a  variational  analysis  of  the  vibrations  of  the  piezoelectric 
thin  film  on  semiconductor  composite  resonator  it  was  shown  that  the 
quality  factor  (Q)  due  to  radiation  into  the  semiconductor  wafer  is  a 
sharply  varying  function  of  the  ratio  of  the  thickness  of  the  resonator 
to  that  of  the  film  if  trapping  is  not  present,  but  not  if  trapping 
is  present.  The  treatment  is  very  cumbersome  to  use  and  is  restricted 
to  the  case  of  strip  electrodes  and  diaphragms.  In  this  work  a 
perturbation  procedure  for  the  calculation  of  the  Q  due  to  radiation 
into  the  wafer  is  presented,  which  is  considerably  less  cumbersome  to 
use  than  the  earlier  treatment  and  is  not  restricted  to  the  case  of 
strip  electrodes  and  diaphragms.  The  resonant  mode  of  interest  is 
determined  from  an  equation  for  transversely  varying  thickness  modes 
in  composite  resonators  and  simple  approximate  but  accurate  conditions 
at  the  edges  of  the  diaphragm,  from  which  the  radiation  into  the  wafer 
is  obtained  using  a  variational  approximation  procedure.  The  resonant 
mode  and  resulting  radiation  field  are  employed  in  a  perturbation 
integral  to  calculate  the  Q.  For  the  case  of  strip  electrodes  and 
diaphragms  the  calculated  results  are  shown  to  be  in  good  agreement 
with  the  earlier  more  cumbersome  calculations.  In  addition  the 
perturbation  calculations  are  performed  for  the  case  of  rectangular 
electrodes  and  diaphragms. 


1.  Introduction 


The  composite  resonator  consists  of  a  uniform  thin  layer  etched 
in  a  small  well-defined  region  of  a  semiconductor  wafer  to  form  a 
diaphragm,  upon  which  is  deposited  a  thin  piezoelectric  film  along  with 
the  electrodes  to  form  a  resonant  region  directly  on  the  wafer.  In  a 
recent  analysis*  of  the  composite  resonator  vibrating  in  an  essentially 
thickness-extensional  mode  it  was  shown  that  the  quality  factor  (Q)  due 
to  radiation  into  the  bulk  semiconductor  is  a  very  sharply  varying 
function  of  the  ratio  of  the  thickness  of  the  diaphragm  to  that  of  the 
film  if  trapping  is  not  present,  but  not  if  trapping  is  present.  A 
brief  explanation  of  the  meaning  of  the  words  "energy  trapping"  is 
given  in  the  Introduction  of  Ref.l.  The  treatment  employs  a  very 
accurate  but  extremely  cumbersome  variational  approximation  technique 
and  is  restricted  to  the  case  of  strip  electrodes  and  diaphragms. 

In  this  work  a  perturbation  procedure  for  the  calculation  of  the  Q 
due  to  radiation  into  the  semiconductor  wafer  is  presented,  which  is 
considerably  less  cumbersome  to  use  than  the  earlier  variational  treat¬ 
ment  and  is  not  restricted  to  the  case  of  strip  electrodes  and  diaphragms 
In  the  treatment  the  resonant  mode  of  interest  is  determined  from  the 

equation  for  transversely  varying  essentially  thickness-extensional 

2 

modes  in  composite  resonators  and  simple  approximate  but  very  accurate 

conditions  at  the  edges  of  the  diaphragm.  This  resonant  mode  is  then 

used  to  determine  the  near  field  radiation  into  the  semi-conductor  wafer 

3 

by  means  of  a  variational  approximation  procedure  .  Of  course,  as  in 
the  earlier  work*,  all  radiating  plate  waves  in  the  wafer  are  included 
in  order  to  achieve  accuracy.  Finally,  the  mode  in  the  composite 


resonator  and  the  radiation  field  in  the  wafer  are  employed  in  a  perturb 


ation  integral  to  calculate  the  Q.  Calculations  are  performed  for  the 


cases  of  rectangular  electrodes  and  diaphragms  and  strip  electrodes  and 


diaphragms  both  when  trapping  is  and  is  not  present.  For  the  case  of 


strip  electrodes  and  diaphragms  the  calculated  results  are  shown  to  be 


in  good  agreement  with  the  earlier  more  cumbersome  direct  calculations  . 


As  in  the  earlier  work  ,  all  calculations  are  performed  for  the  particu¬ 


lar  case  of  an  aluminum-nitride  film  on  gallium-arsenide.  Since  the 


fundamental  essentially  thickness-extensional  mode  will  not  trap  for 


an  aluminum-nitride  film  on  a  gallium-arsenide  diaphragm  in  the  flat 


plate  configuration  and  the  second  thickness-extensional  mode  will 


trap  for  a  sufficiently  thick  aluminum-nitride  film  ,  the  Q  is  calcu¬ 


lated  for  the  second  essentially  thickness-extensional  mode  when 


trapping  is  present  and  for  the  fundamental,  when  trapping  is  not 


present. 


2.  Perturbation  Procedure 


Since  the  coupling  is  small  in  the  piezoelectric  thin  film  and 


the  semiconductor  is  assumed  to  be  nonpiezoelectric,  we  need  consider 


only  the  elastic  portion  of  the  equation  for  the  first  perturbation  of 


the  eigenfrequency,  which  may  be  written  in  the  form 


A  *H/2tu  ,  uu»u)  -  A  , 
A  H  U  M-  n 


(2.1) 


where  and  (u  are  the  unperturbed  and  perturbed  eigenfrequencies , 


respectively,  and 


*V '  I  W"-“iT?j)ds. 


(2.2) 


where  n^  denotes  the  outwardly  directed  unit  normal  to  the  surface  S 
of  the  resonator.  Since,  as  already  noted,  any  piezoelectric  coupling 
is  ignored,  the  constitutive  equation  for  the  stress  tensor  T  is 
given  by 


^ij  ^jkX^hcjX’ 


(2.3) 


where  the  c^^  denote  the  elastic  constants  and  u^  is  the  mechanical 


displacement  vector.  Standard  Cartesian  tensor  notation  is  employed 
along  with  the  summation  convention  for  repeated  tensor  indices  and  the 
comma  convention  for  partial  differentiation  with  respect  to  a  space 


coordinate,  as  in  Ref. 4.  The  normalized  displacement  field  g1?  is 


defined  by 


u^/N  , 


N. 


<f  pu4uj'dV- 


(2.4) 


in  which  u^  denotes  the  purely  real  mechanical  displacement  field 

associated  with  the  mode  of  interest  in  the  composite  resonator  and 

T^.  denotes  the  associated  stress  tensor.  A  cross-section  of  the 
ij 

composite  resonator  attached  to  the  semiconductor  wafer  is  shown  in 
Figure  1.  A  schematic  plan  view  of  the  assumed  composite  resonator  is 
shown  in  Figure  2  along  with  the  assumed  approximate  edge  condition 
(u^ = 0)  for  the  transverely  varying  thickness-extensional  mode, in 
accordance  with  the  explanation  in  the  next  section.  Thus,  for  the 
problem  at  hand  the  surface  S  in  (2.2)  denotes  the  surface  along  which 
the  semiconductor  portion  of  the  composite  resonator  abuts  the  wafer. 
The  calculated  mode  in  the  composite  resonator  shown  in  Figure  2 
results  in  tractions  t^  and  displacements  u^  that  are  applied  to  the 


wafer  along  the  surface  at  which  it  abuts  the  resonator  as  shown  in 


Figure  3.  These  tractions  t^  and  displacements  then  cause  radiation 

fields  in  the  wafer,  which  are  calculated  by  means  of  a  very  accurate 

variational  approximation  procedure  in  Sec. 5.  These  radiation  fields 

in  the  wafer  produce  a  reaction  back  on  the  composite  resonator  along  S. 

The  stresses  T. .  and  displacements  u.  in  (2.2)  denote  this  back  reaction 
ij  J 


field. 


Since  the  radiating  fields  in  the  wafer  are  chosen  to  satisfy 


the  differential  equations  in  and  boundary  conditions  on  the  major 
surfaces  of  the  wafer  exactly,  in  the  purely  elastic  case  all  that 
remains  of  the  appropriate  variational  principle  in  which  all  conditions 
are  natural  conditions  is  given  by^>^ 


I  (tr\iki,6aidS  +  I  “k(W6\.tdS*0’ 


(2.5) 


where  S  and  S  denote  the  portions  of  the  surface  along  the  left  end 
IN  L* 

of  the  wafer  shown  in  Figure  3  on  which  natural-  and  constraint- 
type  conditions  ,  respectively,  are  prescribed.  We  note  that  S^  and  S( 
refer  to  different  portions  of  the  surface  for  different  terms  in  the 
boundary  integrals  depending  on  each  actual  condition  at  a  point.  In 
(2.5)  n^  denotes  the  outwardly  directed  unit  normal  to  the  wafer  and 
we  note  that  for  surfaces  along  which  the  diaphragm  abuts  the  wafer 
the  n  in  (2.2)  and  (2.5)  are  equal  in  magnitude  and  opposite  in  sign. 
In  (2.5)  t ^  and  u^  are  known  from  the  resonant  eigensolution  T*^,  u^ 

A  A 

and  u^  and  1^  are  found  from  the  expressions  for  the  solution  field 

radiating  into  the  wafer,  which  are  obtained  in  Sec. 5.  At  this  point 

*  * 

we  note  that  the  reaction  u  T,  .  from  the  radiating  field  is  complex, 


the  real  part  of  which  yields  a  small  change  in  frequency  when  substi¬ 
tuted  in  (2.2),  which  is  negligible  and  not  of  interest  here,  and  the 
imaginary  part  of  which  yields  the  attenuation  due  to  radiation  into 
the  wafer  from  (2.2)  which  is  sought  here.  From  the  well-known  rela¬ 
tion  m  *  cd^  -  iau1J_/2Q  and  (2.1),  we  obtain 

Q-'l^/V  <2-' 


3.  Transversely  Varying  Thickness-Extensional  Modes 

As  noted  in  Sec. 2,  a  plan  view  for  the  determination  of  the  mode 
of  interest  in  the  composite  resonator  is  shown  in  Figure  2,  in  which 
the  -  denotes  the  electroded  region,  the  S,  the  side  region,  the  T,  the 
top  (and  bottom)  region  and  the  C,  the  corner  regions.  This  notation 
is  essential  for  the  treatment  of  the  trapped  modes,  but  is  not  needed 
for  the  untrapped  modes  because  in  the  untrapped  case  the  edges  of  the 
electrodes  are  relatively  unimportant.  It  has  been  shown  in  Sec.V  of 
Ref. 2  that  the  homogeneous  equation  governing  the  essentially  thickness- 
extensional  modes  may  be  written  in  the  form 


..  /3  fn  .  a  fn\  -f  ~2  ,n  .  f._.n  . 

Mn {—2  +~TJ  -  +  "  f  -°- 


where  n  denotes  the  order  of  the  pure  thickness  mode  in  the  composite 
resonator,  both  the  superscript  and  subscript  f  denote  the  film  and 
Tffn=T1fn  in  an  electroded  region  while  Tj^  *  in  an  unelectroded 


fn  fn 


region,  which  are  defined  in  Eqs.(4.48)  and  (4.46),  respectively,  of 
Ref. 2  with  the  aid  of  the  appropriate  root  of  (3.22)  of  Ref. 2,  which 
gives  T]^n  and  with  (3.20)  T|°n,  and  c^  is  defined  in  Eq.(3.12)  of  Ref.  2. 


6. 


The  expression  for  the  coefficient  M  ,  which  is  very  important  because 
its  sign  indicates  whether  trapping  is  or  is  not  present  in  the  flat 
composite  plate  configuration,  is  given  in  Eq.(4.44)  along  with  (4.41), 
(4.42),  (4.37),  (4.38)  and  (4.33)  of  Ref . 2  and  clearly  is  much  too 

g 

cumbersome  to  present  here  .  It  has  also  been  shown  in  Sec.V  of  Ref . 2 


that  the  dominant  u^-displacement  field  accompanying  the  mode  is  given 


tn  n  .  .  sr 

3u3  f  (xl’x2’C)'  u3 


sn_n.  . 

u  f  (Xl,x„t), 


(3.2) 


where  the  superscripts  f  and  s  denote  the  film  and  semiconductor 


diaphragm,  respectively,  and  qu^  represents  the  thickness  dependence 
only  and  is  given  by 


and  the  amplitudes  are  given  in  Eqs.(5,20)  of  Ref . 2.  Furthermore,  it 
has  also  been  shown  in  Secs.  IV  and  V  of  Ref. 2  that  at  an  interface 
between  electroded  and  unelectroded  regions  we  have  the  continuity 
of  the  two  quantities 

fn,  df°/dn,  (3.4) 

where  d/dn  denotes  the  normal  derivative  (here  either  d/dx^  or 

d/dx^.  However,  the  boundary  conditions  along  the  edge  of  the 

assumed  resonator  shown  in  Fig. 2  that  are  consistent  with  the 

approximation  made  in  obtaining  Eq.(3.1)  must  be  explained. 

It  has  been  shown  that  to  lowest  order  the  pertinent  consti- 

2 

tutive  equations  for  in-plane  tractions  take  the  form 


(3.5) 


T11  =  C13U3,3  »  T22  =  C13U3,3  ’  T13  “  C44(u3, 1  +  U1, 3)  » 


T23  "  C44(u3,2+u2,3)  ’  T12  “  C66(ul,  2  +u2, 1}  =  °  » 


for  both  the  piezoelectric  thin  film  and  the  diaphragm  essentially 
because  u^  and  u2  are  an  order  of  magnitude  smaller  than  u^  for  small 
wavenumbers  along  the  plate  and  each  differentiation  with  respect  to 
or  x2  reduces  the  order  of  magnitude  by  one.  On  an  edge  normal 

9 

to  x^  the  uniqueness  theorem  reveals  that  we  must  specify  one  term 
of  each  of  the  three  products 


T  u  T  u  T  u 
11  1  '  12  2’  13  3  - 


(3.6) 


However,  by  virtue  of  the  aforementioned  ordering  T^2u2  is  two  orders 
of  magnitude  smaller  than  each  of  the  other  two  terms  and,  hence,  may 
be  neglected,  as  already  indicated  in  (3.5)^.  Consequently,  on  a 
free  edge  normal  to  x^  we  should  satisfy 


0  T 
»  13 


(3.7) 


and  on  a  fixed  edge  normal  to  we  should  satisfy 


u  =  0  ,  u  =  0  . 


(3.8) 


Similar  considerations  for  an  edge  normal  to  x2  reveal  that  on  a  free 


edge  we  should  satisfy 


T  =0  T  -  0 


(3.9) 


and  on  a  fixed  edge  we  should  satisfy 


u2  -  0  ,  u3  =  0 


(3.10) 


Since  in  the  approximation  made  in  obtaining  Eq.(3.1)  governing 
the  essentially  thickness -ex tens ional  modes,  we  have  eliminated  all 
waves  except  the  important  one,  we  can  satisfy  only  one  of  the  two 
conditions  in  each  of  (3.7)-  (3.10).  However,  since  for  small  wave- 
numbers  along  the  plate  is  large  while  and  u^  are  small,  from 
this  fact  and  (3.5)  it  is  clear  that  one  of  the  two  conditions  in  each 
of  (3.7)-  (3.10)  is  large  and  the  other  is  small.  Consequently,  we 
take  the  small  equation  in  each  of  (3.7)-  (3.10)  to  be  satisfied 
approximately  and  require  the  solution  to  satisfy  the  large  equation 
in  each  of  (3.7)-  (3.10)  only.  Since  in  each  instance  the  large 
term  in  each  of  (3.7)  -  (3.10)  is  either  u^  or  u^  we  have  shown  that 
for  either  completely  free  or  completely  fixed  conditions  along  the 
edges  of  the  diaphragm  shown  in  Fig. 2,  the  appropriate  condition  is 

f“  -  0  ,  (3.11) 

on  account  of  (3.2)  and  (3.3).  From  Fig.l  it  is  clear  that  the  edge 
of  the  film  is  free  while  the  edge  of  the  diaphragm  is  essentially 
fixed,  and  we  have  shown  that  in  this  approximation  (3.11)  is  the 
appropriate  condition  for  either  case.  This  means  that  within  this 
approximation  both  u^  and  either  or  T^2  vanish  simultaneously 
along  the  edge  of  the  resonator. 

Since  only  the  semiconductor  diaphragm  abuts  the  wafer  as  shown 
in  Fig.l  only  the  variables  in  the  semiconductor  are  relevant  here. 
However,  in  order  to  calculate  the  tractions  t  and  displacements  u 

At  At 

that  the  mode  in  the  composite  resonator  exerts  on  the  wafer  as  shown 
in  Fig. 3,  we  need  the  associated  displacement  fields  u^  and  u^  in  the 


9. 


semiconductor  diaphragm  in  addition  to  u^. 

In  the  untrapped  case  for  the  fundamental  mode  for  any  n,  f11 
is  given  by 

..n  —  —  iuut 

f  -  cos  §x^  cos  vx2  e  ,  (3.12) 

the  substitution  of  which  in  (3.1)  yields 

f  2  _f  -2  —2  —2 

pV-c^+M^+v  ),  0.13) 

in  which  for  convenience  we  have  assumed  that  the  upper  electrode 
covers  the  entire  film  since  the  edges  of  the  electrode  are  unimportant 
because  they  do  not  cause  exponential  decay  in  the  unelec troded  region 
in  the  untrapped  case.  Substituting  from  (3.12)  into  the  edge  condi¬ 
tions  (3.11)  for  the  mode  of  interest,  we  obtain 

5  *  rr/2d  ,  v  *  rr/2w  ,  (3.14) 

with  which  the  unperturbed  resonant  frequency  cu  may  be  obtained  from 
(3.13).  It  has  been  shown  in  Sec.V  of  Ref . 2  that  for  (3.12),  we  have 

u®  *  tS(x3)  sin  gx1  cos  vx2  eitut  , 

uS2  *  ts(x3)  cos  ^  sin  vx2  eiu,t  ,  (3.15) 

where 

tS  -  A*1  sin  T|°x3+A®2  sin  cos  1)^ 

+  B®2  sin  hS71°x3  ,  C2  *  §2  +  ^  > 


(3.16) 


10 


and  h  and  the  amplitudes  are  given  in  (4. 30 )  and  (4. 26),  respectively, 
of  Ref. 2. 

In  the  trapped  case  we  have  Eq.(3.1)  with  T)  ■  Tj,  in  the 

tn  rn 

a 

electroded  region  denoted  _  and  T|  ■  *-n  the  unelectroded  regions 

denoted  S  and  T  in  Fig. 2  and  as  usual  ignoring  (3.1)  in  the  relatively 
unimportant  corner  regions  denoted  C.  In  addition  the  solution  must 
satisfy  the  continuity  conditions  in  (3.4)  between  the  —  and  S  and  T 
regions  and  between  the  S  and  C  and  T  and  C  regions  and  the  conditions 
(3.11)  along  the  edge  of  the  diaphragm,  i.e.,  at  x^  =*±d  and  x2=±w. 
Clearly  then  for  the  symmetric  modes  in  x^  and  x^  the  expressions  for 
fn  in  each  of  the  respective  regions  shown  in  Fig. 2  are  given  by 


__  _  _  o  -^(x^i)  §(x  -jO 

t  *  E  cos  §x^  cos  vx 2  ,  r  *  (E_e  +  E+e  )  cos  vx2  , 


,nT  .  T  ~v<x2”b)  T  v(x„-b)  - 

f  *  (E_e  +  E+e  2  )  cos  §x^  , 

nC  C  -§(x1-^)-v(x2-b)  c  -§(x  -X)*v(x  -b) 
f  -  Exe  +E2e 

c  §(x1-i)-v(x2-b)  c  §(x1-£)+v(x2-b) 

+  E_e  +  E.  e  ’ 

3  4 


the  substitution  of  which  in  (3.1)  yields  (3.13)  along  with 


(3.17) 


f  2  -f  *2  M  .  2  -2, 
p  *  ’  c331'fn  '  V5  ‘v  >■ 

£  2  _£  -2  ,  2  -2. 

p  "  '  C33V  '  Mn(v  ‘  5  >  • 


(3.18) 


Substituting  from  (3.17)  into  (3.4)  at  x^i  and  x2  -  b  and  into  (3.11) 
at  x^  ®  d  and  x2  =  w,  we  obtain 


11. 


§  tan  - 1 =  c 


7..,  (l+e-2^) 


1  -  e 


-2§(d-i)  ’ 


v  tan  vi =  v 


fl+e-2v(«-t,)> 
-2v(w-b)  ’ 


(3.19) 


1  -  e 


where 


ES 


E  cos  Zl  T 

’  E± 


E  cos  vb 


±2?(d-2) 


1  -  e  ^  “  1  -  e 

,C  E  cos  EJj  cos  vb 


±2v(w-b)  ’ 


1  (i.e'2'(d'£))(l.e‘2v(w‘b))  ’ 


^..e-ZvCw-b)^  EC  _e-25(d.Z)ECf 


C  =  - 2|(d-£)-2v(w-b)  C 

6  h  • 


(3.20) 


Now,  as  in  Ref. 2,  from  (3.13),  (3.18)  and  Eqs.(3.21),  (3.23), 
(3.31)  and  (3.32)  of  Ref .2,  we  obtain 
.f 


r2c"  n  9  .  --,1/2  2cJ.  -  An  9-.l/ 2 

33  ~n  —2  I  33  ,„oN 2  *n  —2 

5  *  LIT"  (V  A  ■  5  J  ,  V  -  [—  Olf)  A  -  v  J  ,  (3.21) 


where 


ln  =  .  (P0n  +R")/G°n 


(3.22) 


and  R",  G°n  and  P°n  are  defined  in  Eqs.(3.19),  (3.24)  and  (3.33), 


respectively,  of  Ref. 2.  Equations  (3.19),  with  (3.21)  and  (3.22) 
constitute  two  independent  transcendental  equations  for  §  and  v  for 
a  given  £,  b,  d  and  w,  which  may  readily  be  solved  for  the  fundamental 
mode  of  interest.  The  eigenfrequencies  for  that  mode  may  then  be 
determined  from  (3.13). 


■M 


•3 


<3 


•'.1 


Nn 


& 


. •r.  ^ 


In  the  trapped  case  since  we  have  growing  and  decaying  solu¬ 


tions  in  the  S  and  T  regions,  from  Secs. Ill  and  V  of  Ref. 2  it  is 
clear  that 


sS 


sT 


C 


s  t"(x3)[e> 


■S(x.-l) 


?(*,-*) 


-  ] 


cos  vx , 


ji  tS(V[E- 


T  v(x2~b)  T  v(x2*b) 
e  -  E+e 


]  cos  ?x:  , 


(3. 23) 


where 


(O 


S.  2  -2 


v  -  ?  ,  (C  ) 


T.  2  —2 


(3.24) 


,  sS  sT 

and  u^  is  for  the  side  region  and  u2  is  for  the  top  region.  We 

further  observe  that  since  all  wavenumbers  along  the  resonator  are 
g  s 

small,  u2  in  the  S  region  and  u^  in  the  T  region  are  negligible  for 

radiation  in  the  x^-  and  x2-directions,  respectively,  in  the  wafer, 

as  discussed  more  completely  in  Sec. 4.  Moreover,  it  is  clear  that 

sC  sC 

for  the  corner  regions  u^  and  u2  may  be  obtained  from  (3.23)  simply 

by  replacing  the  respective  trigonometric  functions  in  each  term  by 

C 

the  exponential  function  in  brackets  in  the  other  term  and  '  is 
given  by 


(CG)2  =-  (§2+v2)  . 


(3.25) 


4.  Variable-Crested  Waves  in  Wafer 

In  this  section  we  obtain  the  solution  functions  for  the  near 
field  waves  with  slowly  varying  crests  and  the  associated  dispersion 
relations.  The  displacement  equations  of  motion  and  the  linear  elastic 
constitutive  equations  for  the  nonpiezoelectric  gallium-arsenide 
semiconductor  wafer  with  x_  along  a  cube  axis  take  the  form^ 


At 


13. 


cll“l,U  +  <C13+C44,“2,12  +  (<:n  +c44)'*3,13  +  c44"l,  22  +C44S1,33  ' 
c44“2, 11  +  <c44  +C13)G1,12  +  Clia2,22  +  (c13  +  C44)S3,  23  +C44G2,33  '  =52 
c44“3,U  +  <C44+c13);:i,31+';4433,22  +  (c44+c13)52,23  +C3323,33*=53 


(4.1) 


T11  =  C11Q1,1+C13U2,2+C13U3,3  ’  T22  “  C13U1, 1  +  C11U2,  2  +  C13U3, 3  » 


T33a‘C13Ul,l+C13U2,2+C33U3,3’ 


T23  =  C44(G3,2+G2,3)’  T12  *  c44(fil, 2  +  *2, 1}  ’  T31  “  C44(u3, 1  +  ’ ul,  3  ) 


(4.2) 


where  the  carets  are  used  to  denote  the  wafer  as  distinct  from  the 


diaphragm. 


From  either  Figs.l  or  3  we  see  that  the  boundary  conditions 


on  the  major  surfaces  of  the  semiconducting  plate  may  be  written  in 


the  form 


T*  .  =  0  at  x*  0  and  x  =  -  hs  . 
3j  3  3 


(4.3) 


In  considering  waves  radiating  in  the  x^-direction,  we  first  note 


that  since  v  and  v  are  small,  from  (4.1)  and  (4.2)  it  is  clear  that 


T12  is  an  order  of  magnitude  smaller  than  and  u®,  and  u^  is  an 

A  S  A  S  A  C 

order  of  magnitude  smaller  than  u^  and  u^.  Hence,  u^,  the  x .,-di f ferent ial 


equation  and  are  negligible  for  radiation  in  the  x^-direction. 


Under  these  circumstances,  from  (4.1)  we  see  that  the  differ¬ 


ential  equations  that  must  be  satisfied  take  the  reduced  form 


s  r,s  ,  ,  s  ,  s  ,s*s 

CU1,1,11  +  <C13+C44)U3,13+C44U1,33 


S*  s 
0  ux  , 


S  *  S  .  .  S  ,  S  x  “S  S  AS 

C44U3,11  +  <C13+<:44)U1,13+C33U3,33 


S-S 
0  u3  , 


(4.4) 


V.-'.'  .-'.-f.  v. 


h 


14. 


and  from  (4.2)  and  (4.3)  the  boundary  conditions  take  the  reduced  form 


s  -s  .  s  -s 

:13ul,l+c33"3,3’° 

at  =  0 

and  x^ 

=  -hs , 

(4.5) 

s  ,-s  -s  ,  . 

C44(U1,3+U3,1)'° 

at  x^  =  0 

and  x^ 

v  s 

■  -  h  . 

(4.6) 

For  waves  with  slowly  varying  crests  in  the  x^direction  and 
propagating  in  the  +x1"direction  in  the  near  field  consider 


*s  -  “S  iTVS  -S  “iT1SX3  i§<Xrd)  imt 

ua  =.  cos  vx2(Aae  3+B*e  3)  ,  e  e1UJt, 

which  satisfies  (4.4)  provided 


a  *  1,3  ,  (4.7) 


where 


s  ts  ,  s  ?s 


ailAl  +al3A3  '°  ’  al3Al  +CT33A3  ”  °  ’ 

0,  •M+a'&.O, 


S  S  *  2  e  ^  2  c  2 

ali  -cn5  -cJ41C  +  p  »  . 


,  s  s  .  %  - 
(C13  c44)5  71  s  » 


_s  s  %2  s  -2  ,  s  2 

33  =  '  C44'=  "  C33%  p  ^  • 


(4.8) 

(4.9) 


(4.10) 


Each  of  the  two  systems  of  linear  homogeneous  algebraic  equations  in 

two  amplitudes  yields  nontrivial  solutions  when  the  determinant  of  the 

coefficients  of  the  amplitudes  vanishes.  Both  determinants  are 

-2  *2  2 

identical  and  each  is  quadratic  in  §  ,  "l  and  w  .  Thus,  for  a  given  • 

A  »(i)  »(2). 

and  ju  each  determinant  yields  two  independent  7]  (T|g  ,  "r|g  )  and  each 
7]^^  yields  amplitude  ratios  from  each  system  of  linear  algebraic 
equations,  i.e.,  (4.8)  and  (4.9).  Hence,  there  are  four  amplitude 

ratios,  which  we  denote  by 

•si  si-si  jsi  si-si 

A  -  u  A  B,  =  v  B 


(4.11) 


15. 


As  a  solution  of  the  boundary  conditions  (4.5)  and  (4.6),  we 


u^  =  cos  vx2 


J;  iTl  x,  -iTl  x.  i§(x.-d) 

V  ^sm  'sm  3  ,  £sm  stn  3.  5  1  i< 

)  (A3  e  +B3  e  )  e  e 


u^  =  cos  vx2 


I  <4 


iTl  x„  -ill  x.  i|(x,-d)  . 

sm»sn  'sm  3  ,  Asm£sm  sm  3.  1  lcut 

A^  e  +  v  e  )e  e 


(4.12) 


Substituting  from  (4.12)  into  the  boundary  conditions  (4.5)  and  (4.6) 


we  obtain 
2 


Z”?sn.  s  **sn  ,  s  *  .  t  ‘sn,  s  c-  sn  s  «,  .1  _ 

_a3  (c13?M.  +  c33^sn  +  B3  (c13^V  “  C33‘1sn)J  =  0  » 


0  ^  A 

-  -  iT|  hs  .  ^  ^  iT]  hs 

V  |Asn  s  pAsn  s  A  .  sn  ,  *sn.  s  JAsn  s  A  .  sn 

L  LA3  (c135“  +c33\n)e  +B3  (c135V  '  c331W)e 


3  '13- 


'33  sn 


t- 

n-sn  -  *5n  , 

A3  ^sn11  +  ?> 


lT|  hs  IT]  hs- 

e  S"  +53n<-\n0  +5>«  S”  ]-0. 


(4.13) 


Equations  (4.13)  constitute  a  system  of  four  linear  homogeneous  alge- 

a  sn  a  sn 

braic  equations  in  A^  and  ,  which  yields  nontrivial  solutions  when 

the  determinant  of  the  coefficients  vanishes.  Calculations  are 

12  a  *sn  Asn  Asn  Asn 

performed  in  the  usual  way  ’  and  yield  the  |i  ,  v  ,  A  and  B 

A 

along  with  the  dispersion  relation  <ju=u;(§). 

The  functions  in  (4,12)  are  for  either  the  untrapped  case  or 
the  electroded  region  for  the  trapped  case  and  since  for  slowly  varying 
crests  only  the  slowly  varying  terms,  i.e.,  the  cos  vx?  in  (4,12), 
which  have  no  influence  on  the  radiating  waves  or  the  dispersion 
relation,  change,  it  is  not  purposeful  to  write  the  solution  for  any 
other  regions  for  the  trapped  case.  From  the  solution  the  lowest  17 


real  dispersion  curves  for  the  bulk  gallium-arsenide  plate  have  been 


calculated  and  are  shown  in  dimensionless  form  in  Fig. 4,  which  is 
identical  with  Fig. 6  of  Ref.l  since  for  slowly  varying  crests,  the 

A 

crests  do  not  influence  the  dispersion  relation  ud  =  'ju(I).  For  the 
fundamental  essentially  thickness-extensional  resonance  of  the 
composite  resonator  consisting  of  a  7  pm  thick  aluminum-nitride  film 
on  a  14  um  thick  gallium-arsenide  diaphragm,  which  is  around  132  MHz, 
this  number  of  dispersion  curves  is  for  a  gallium-arsenide  wafer 
thickness  of  about  5  mils.  In  this  work  we  perform  calculations  for 
wafer  thicknesses  up  to  8  mils  for  which  there  are  30  real  dispersion 
curves  for  a  frequency  of  132  MHz.  However,  we  do  not  bother  to  show 
the  figure  for  more  than  17. 

Similarly,  since  statements  equivalent  to  the  foregoing  hold 
for  waves  with  slowly  varying  crests  in  the  x^-direction  and  propagating 

in  the  X2~direction  in  the  near  field,  the  differential  equations  and 
boundary  conditions  that  must  be  satisfied  may  be  obtained  from 
(4.4)  -  (4.6)  simply  by  replacing  all  indices  1  by  2  in  (4.4)  -  (4.6). 
Under  these  circumstances  the  near  field  radiating  solution  takes  the 


i;  ill  x.  -it)  x_  iv(x„-w)  . 

As  7  v  ,?sm  sra  o  .gsm  'sm  3x  2  i< 

u^  =  cos  Sx1  )  (A ^  e  +B3  e  )e  e 

m=l 


-s 

u„  =  cos 


—  V  *sm*sm  ^\mX3  Asm*sm  ^smX3<.  iV^X2  i^t 
cx.  )  (ll  A_  e  +v  B„  e  )e  @ 


in  place  of  the  form  shown  in  (4.12).  From  this  solution  we  obtain  the 


propagating  dispersion  curves  for  near  field  radiating  waves  in  the 
x^-direction.  The  dispersion  curves  are  identical  with  those  in  Fig. 4. 


5 .  Variational  Analysis  of  Radiation  into  Wafer 

In  this  section  we  determine  the  waves  radiating  into  the  semi¬ 
conductor  wafer  due  to  the  tractions  and  displacements  resulting  from 
the  mode  in  the  composite  resonator  by  means  of  a  variational  approxi¬ 
mation  procedure,  as  indicated  in  Sec. 2.  Since  the  near  radiation 
fields  emanate  from  the  edges  of  the  diaphragm,  which  are  normal  to  x^ 
and  x^,  respectively,  we  may  determine  the  near  fields  radiating  in 
the  x^-  and  X2*-directions  separately.  Accordingly,  we  first  consider 
radiation  in  the  x^direction.  From  Fig.  3  we  see  that  the  edge  of  the 
wafer  below  the  diaphragm  is  traction  free,  i.e.,  T^=0.  In  Sec. 3  we 


have  shown  that  for  the  mode  in  the  resonator  vanishes  to  the 


order  of  approximation  throughout  and  u^  and  both  vanish 
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simultaneously  along  the  edge  of  the  resonator  ,  i.e.,  where  the 
diaphragm  abuts  the  wafer,  while  both  and  u^  exist  along  the  same 
edge.  Furthermore,  in  Sec. 4  it  is  shown  that  since  the  transverse 


mode  shape  is  slowly  varying  in  the  wafer,  for  radiating  waves  propa¬ 


ll 


gating  in  the  x^-direction  in  the  wafer  u^  is  negligible  .  As  a 


consequence  of  the  foregoing,  for  propagation  in  the  +x^-direction 


in  the  untrapped  case  Eq.(2.5)  takes  the  form 


18. 


I  dx2[J  "n6Sl+f13653)l  *,-ddx3+  J  ^U66! 


-w  -s 

-h 


+  (-T®3+T*3)6UjJ|  d  dx3  +  J  [-(SJ-^STJ, 

-1-  q 

-hS 

-  G3«S13]I, ,-d  d*3]  * 


(5.1) 


where,  as  already  indicated,  the  forcing  terms  resulting  from  the  mode 

s  s 

in  the  composite  resonator  are  and  g^-  In  the  trapped  case,  the 

integral  from  -w  to  w  is  replaced  by  integrals  over  the  S  and  C  regions, 

respectively,  i.e.,  from  -w  to  -b  and  -b  to  b  and  b  to  w,  in  accordance 

S  S 

with  the  solution  in  (3.17)  and  the  expressions  for  Uj  and  u 2  given 
and  discussed  at  the  end  of  Sec. 3.  For  radiation  in  the  x2~direction 
we  simply  interchange  subscripts  1  and  2  and  replace  w  by  d  in  (5.1) 
for  the  untrapped  case  and  w  by  d  and  b  by  L  for  the  trapped  case. 

For  radiation  in  the  x^-direction  we  now  expand  the  solution 
in  the  wafer  as  a  sum  of  waves  with  slowly  varying  crests  in  the 
x2“direction,  which  were  discussed  in  Sec. 4  and  are  given  in  (4.12). 


£(a)*s (a) 
KS  Ua 


a-  1,3, 


(5.2) 


-  '  .*  V  V  >  ’•> ' 


*  y^'  yv  yv 


fe: 
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where  each  of  the  uS ^  are  of  the  form  given  in  (4.12)  and  N  denotes 

cL 

A 

the  number  of  branches  of  the  dispersion  curves  required  for  the  X 

,  ,  -s (a) 

propagating  plate  waves  for  a  given  wafer  tnickness.  Since  the 
are  fixed  only  the  are  varied  when  (5.2)  is  substituted  into 

(5.1) .  Accordingly,  substituting  from  (5.2)  into  (5.1),  employing 

(4.2)  and  performing  the  integrations,  we  obtain  an  equation  of  the 

form 


N  N 


I  [I  «,<%+•$]  «sw-». 

3=1  o=l 


(5.3) 


where 


ZP  r  ,*a  *8  so  ,*sm*sm  ,  -sm.  jsn 

I  [*<1C.+15n,-i.<Vg  +  T5=2  A3B 


m=l  n=l 


~a  2,3  ,o  .-sm-sn  ,  *sn,  2sm 


-o  20  .o  -sm^sn  ,  -sn.qsn  .“or  2,\°  "sm^sn  Asn.»sm 

+  S'tm-\n>-hs(Tlav3  +T5»)B39’g(T|s».-1')-ns(^  aU  +  053M3c 

‘a  *9  N o  ,*sin*sn  -sm.  «sn  _/  201  ,2.9  \°  ,.Asro*sn  ^  Asnx£sm 


»of  *9  .  o  ,''sra*sn  ,  ~sm.Asn  .  20  ,22  .o  ."sro-sn  -sn,"ssi 
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and  is  the  normalization  factor,  which  is  defined  in  (2.4)  and  for 
the  mode  in  the  composite  resonator  is  given  by 
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In  (5.4)  and  (5.5)  we  have  employed  the  definitions 
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The  foregoing  is  for  the  untrapped  case.  For  the  trapped  case  is 

S 

unchanged  and  may  be  obtained  from  the  expressions  in  (5.5)  simply 
by  replacing  §  sin  §d  by  §  cos  §2/sinh  §(d-£)  and  £  by  £s  and  N  takes 
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in  place  of  (5.6).  Since  the  variations  in  (5.3)  are  arbitrary,  we  obtain 

*  (or) 

the  inhomogeneous  linear  algebraic  equations  for  the  K  in  the  form 
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In  a  similar  way  for  radiation  in  the  X2~direction,  we  obtain 

A 

1  6-l> 

Y=1 


where  Ays  is  given  in  (5.4)  and  may  be  obtained  from  the  express 
in  (5.5)  for  the  untrapped  case  and  the  revised  expression  for  the 
trapped  case  simply  by  replacing  §  by  v,  d  by  w  and  i  by  b. 

A 

Equations  (5.12)  and  (5.13)  each  constitute  N  inhomogeneous  linear 

A  A  /y  ) 

algebraic  equations  for  the  N  unknowns  Kg  from  (5.12)  and  the 
N  unknowns  from  (5.13),  respectively.  When  the  Kg  and  the  K 

have  been  determined,  the  near  field  radiating  solution  is  known. 


6.  Oualitv  Factor  Resulting  from  Radiation 

In  this  section  we  calculate  the  Q  due  to  radiation  into  the 
semiconductor  wafer  from  the  analyses  presented  in  the  previous 
sections.  In  accordance  with  Secs.3-5  it  is  clear  that  in  the 
untrapped  case  the  perturbation  integral  takes  the  form 
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while  in  the  trapped  case  the  integral  from  -w  to  w  is  replaced  by 
integrals  from  -w  to  -b  and  -b  to  b  and  b  to  w  and  the  integral  from 
-d  to  d  is  replaced  by  integrals  from  -d  to  -l  and  -i  to  JL  and  l  to  d 

5 

in  accordance  with  the  solution  in  (3.17)  and  the  expressions  for  u^ 
and  u^  given  and  discussed  at  the  end  of  Sec. 3.  Substituting  the 


imaginary  parts  of  the  solutions  determined  in  Secs. 3  and  5 
into  (6.1)  for  the  untrapped  case  and  into  the  equivalent  integral 
for  the  trapped  case  and  performing  the  integrations,  we 
obtain 


N  „(<*) 

^  -  r  [« ?<»  +  )  I  +  -=r)] 

^  a=l  c 

—  N  u 

+v  3i„  v»(d  l  «fa)(H(.)+  -i-)] 

5  o-l  c 


in  the  untrapped  case  and 
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in  the  trapped  case,  where 
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Either  Eq.(6.2)  or  (6.  3)  may  now  be  used  to  calculate  the  Q  from  (2.6) 
Of  course,  as  in  the  earlier  work1  all  radiating  plate  waves  in  the 
thick  region  of  the  gallium-arsenide  are  included  to  achieve  accuracy. 
Since  at  a  given  frequency  the  number  of  radiating  waves  in  a  plate 
goes  up  significantly  with  thickness,  when  trapping  is  not  present  we 
have  considered  gallium-arsenide  wafers  no  thicker  than  8  mils  at  a 
frequency  around  132  MHz,  for  which  there  are  30  radiating  plate 
waves.  Results  are  presented  for  wafer  thicknesses  ranging  from 
4  rails  to  8  mils  because  this  is  considered  to  be  within  the  practical 
range.  For  the  untrapped  case  the  results  presented  are  for  a  film 
thickness  of  7  microns  and  a  diaphragm  thickness  of  14  microns.  Both 
strip  diaphragms  600  microns  wide  and  square  diaphragms  with  lateral 
dimensions  of  600  microns  x  600  microns  were  considered.  The  strips 
are  treated  for  comparison  with  the  earlier  work1.  The  results  can 
readily  be  obtained  for  the  strip  case  from  the  analysis  for  the 


rectangular  case  simply  by  allowing  w  and  b  to  get  very  large  compared 
with  JL  and  d,  for  which  v  -  0  and  the  C  and  T  regions  are  eliminated 
as  is  radiation  in  the  X2~direction.  For  the  trapped  case  the  results 
presented  are  for  a  film  thickness  of  12  microns  and  a  diaphragm 
thickness  of  14  microns  because  the  trapping  is  considerably  better 
for  this  combination  of  thicknesses  and  trapping  does  not  exist  for 
the  other  at  the  second  thickness  mode^.  For  this  combination  of 
thicknesses  the  second  thickness  mode  is  around  250  MHz.  In  the 
trapped  case  both  strip  electrodes  500  microns  wide  and  square 
electrodes  with  lateral  dimensions  of  500  microns  x  500  microns  were 
considered.  The  lateral  dimensions  of  both  the  strip  and  square 
diaphragms  was  varied  and  the  wafer  thickness  is  6  mils. 

In  the  absence  of  trapping  for  the  strip  case  the  values  of  Q 
calculated  from  Eq.(2.6)  of  this  work  are  plotted  as  the  solid  curve 
in  Fig. 5,  In  which  the  dotted  curve  from  Fig. 8  of  Ref.l  is  also 
plotted  for  purposes  of  comparison.  It  can  be  seen  from  the  figure 
that  the  agreement  is  quite  good.  However,  although  the  highest  Q's 
calculated  in  this  work  are  very  nearly  the  same  as  those  obtained 
in  the  earlier  more  cumbersome  direct  calculation*,  the  lowest  Q's 
calculated  by  means  of  the  perturbation  procedure  tend  to  be  nearly 
an  order  of  magnitude  higher  than  those  calculated  by  the  earlier 
direct  procedure*.  We  are  not  absolutely  sure  of  the  reason  for  this 
discrepancy,  but  there  are  two  possibilities.  The  perturbation 
procedure  might  be  tending  to  lose  its  accuracy  for  low  Q  because  of 
the  increased  radiation  or  the  resonant  frequency  might  not  have  been 
sufficiently  precisely  determined  by  means  of  the  earlier  direct 
procedure*  for  the  accurate  determination  of  the  lowest  Q  values. 
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However,  Fig. 5  reveals  that  the  highest  Q's  calculated  by  the  perturba 
tion  procedure  are  consistently  slightly  higher  than  those  calculated 
using  the  earlier  direct  procedure,  which  tends  to  support  the  second 
possibility.  The  figure  also  shows  that  the  location  of  the  peaks  and 
valleys  of  Q  with  wafer  thickness  determined  by  means  of  the  perturba¬ 
tion  procedure  is  in  quite  good  agreement  with  those  obtained  from 
the  earlier  direct  calculation  .  Also, in  the  absence  of  trapping 
the  Q's  calculated  from  Eq.(2.6)  for  the  square  diaphragm  are  plotted 
in  Fig. 6.  It  can  be  seen  from  the  figure  that  the  peaks  and  valleys 
are  in  essentially  the  same  positions  as  in  the  strip  case,  but  that 
the  Q's  are  considerably  lower,  roughly  between  1/2  to  2/3  of  the 

values  in  the  strip  case.  This  is  as  expected  because  of  the  radia¬ 
tion  in  two  orthogonal  directions  for  rectangular  diaphragms.  When 
trapping  is  present  the  Q's  calculated  from  Eq.(2.6)  for  both  the 
strip  and  square  case  are  plotted  in  Fig.  7  as  a  function  of  (d  -  Z) , 
i.e.,  the  distance  from  the  edge  of  the  electrode  to  the  edge  of  the 
diaphragm.  It  can  be  seen  from  the  figure  that,  as  expected,  Q 
increases  very  rapidly  with  (d  -  Z)  and  for  the  same  value  of  (d  -  Z) 
the  Q  is  about  twice  as  large  in  the  strip  case  as  in  the  rectangular 
case.  If  a  film  thickness  of,  say,  8  microns  had  been  employed,  the 
required  (d  -  Z)  for  good  Q  would  be  much  larger  because  although 
trapping  is  present  for  that  ratio  of  thicknesses,  the  dispersion 

curve  for  the  trapped  mode,  i.e.,  the  value  of  M  ,  indicates  that  the 

7  n 

spatial  decay  rate  at  resonance  is  much  slower  than  for  the  ratio  of 
thicknesses  considered  here. 
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FIGURE  CAPTIONS 


Figure  I  Cross-Section  of  a  Composite  Resonator 

Figure  2  Plan  View  of  Model  for  Analyses  of  Composite  Resonator 

Mode  Shapes 

Figure  3  Cross-Section  for  Variational  Analysis  of  Radiation 

into  Wafer 

Figure  4  Dispersion  Curves  for  the  Gallium -Arsenide  Wafer 

with  Q  the  Dimensionless  Frequency  Normalized  with 
Respect  to  the  First  Thickness-Shear  Frequency 

Figure  5  Quality  Factor  Versus  Wafer  Thickness  when  Trapping 

is  not  Present  for  the  Strip  Composite  Resonator. 

The  dotted  curves  are  from  Fig. 8  of  Ref.l. 

Figure  6  Quality  Factor  Versus  Wafer  Thickness  when  Trapping 

is  not  Present  for  the  Rectangular  Composite  Resonator 

Figure  7  Quality  Factor  Versus  Distance  from  Edge  of  Electrode 

to  Edge  of  Diaphragm  when  Trapping  is  Present 


